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We  have  re-reduced  the  Southern  Polar  Zone  of  the  Yale  Photographic  Star 
Catalogue  (-70°  to  -90°  declination)  using  an  overlapping  plate  technique.  This 
region  was  photographed  in  1955-56  with  64  overlapping  plates  and  the  reduction 
of  these  plates  (i.e.  the  determination  of  the  equatorial  coordinates  of  the  star 
images  on  the  plates)  was  completed  in  I97L  Because  of  the  scarcity  of 
observations  from  the  Southern  Hemisphere  prior  to  the  1970's,  this  data  set 
is  particularly  valuable. 

The  plates  where  photographed  so  that  they  overlapped  50%  in  right  ascension 
and  50%  in  declination.  Thus  every  region  of  the  sky  was  covered  by  at  least 
2  plates  and  sometimes  as  many  as  20  different  plates.   Originally  these  plates 


were  individually  reduced.  Thus  after  the  reduction,  some  of  the  stars  had  many 
as  20  different  position  estimates  (i.e.  right  ascension  and  declination).  Clearly 
a  star  occupies  one  position  at  one  time,  but  this  can  not  be  obtained  without 
mathematically  enforcing  this  constraint  on  a  multiple  plate  reduction  process. 
Instead,  the  separately  computed  values  were  averaged  to  give  the  "best"  (in  the 
mathematical  statistics  sense)  estimate  for  the  stellar  position. 

Taking  advantage  of  the  high  degree  of  overlap  of  these  plates,  we  used  a 
more  powerful  reduction  method  for  overlapping  plates  called  the  overlapping 
plate  technique.  This  method  enforces  the  simple  fact  that,  at  a  given  instant 
of  time,  a  star  is  at  only  one  position.  When  this  constraint  is  enforced,  all  the 
reduction  parameters  on  all  plates  are  obtained  simultaneously  by  solving  one  large 
system  of  equations.  The  application  of  these  improved  reduction  procedures  and 
resources  to  the  original  measurements  have  yielded  star  positions  which  are  a 
20%  improvement  over  those  originally  published.  These  new  estimates  of  star 
positions  are  capable  of  yielding  improved  proper  motions,  which  are  invaluable 
data  for  many  aspects  of  astronomy. 


CHAPTER   1 
INTRODUCTION 

Astrometry,  one  of  the  oldest  branches  of  astronomy,  deals  with  the  determi- 
nation of  locations  and  the  positional  changes  of  celestial  bodies.  The  compilation 
of  star  catalogs  is  the  domain  of  astrometry  devoted  to  the  determination  of  ac- 
curate and  precise  positions  of  celestial  bodies.  One  of  the  most  important  tasks 
of  astrometry  is  to  establish  a  reference  frame  and  maintain  it  by  determining  the 
motions  of  the  reference  stars  which  it  represents.  Usually  these  reference  frames 
are  defined  by  a  star  catalog  (i.e.  FK4  (Fricke  et  al.,  1963),  FK5  (Fricke  et  al. 
1988)  and  soon  the  HIPPARCOS  catalog).  The  compilation  of  star  catalogs  has 
been  a  continuing  processes  over  the  centuries.  One  of  the  first  recorded  star 
catalogs  was  compiled  by  Hipparchus  (128  B.C.)  and  contained  850  stars  (Abell 
et  al.,  1987).  Throughout  the  centuries  astronomers  have  devised  and  improved 
upon  the  methods  to  compute  accurate  and  precise  positions  for  as  many  stars  as 
possible.  Today  the  largest  star  catalog  is  the  Guide  Star  Catalog  (GSC,  Lasker 
et  al.,  1990)  containing  roughly  20  million  objects  with  a  limiting  magnitude  of 
15"^  and  positional  standard  errors  from  the  plate  centers  to  the  edges  which  vary 
from  0".5  to  I'M  in  the  northern  celestial  hemisphere  and  from  T'.O  to  l."6  in  the 
southern  celestial  hemisphere  (Taff  et  al.,  1990).  Containing  fewer  stars,  but  with 
a  much  higher  expected  positional  accuracy,  the  HIPPARCOS  catalog  is  the  first 
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star  catalog  to  contain  observations  solely  made  in  space.  This  catalog  is  being 
produced  by  the  European  Space  Agency  from  observations  taken  by  the  satel- 
lite HIPPARCOS  and  will  contain  roughly  120,000  star  positions  with  limiting 
visual  magnitudes  between  IJ"^  and  8.7"^  depending  upon  the  galactic  latitude. 
The  individual  rms-errors  (root  mean  square-errors)  are  expected  to  be  1.3  milli 
arcseconds  (mas)  per  coordinate  in  position  (Kovalevsky,  1995). 

Star  positions  at  only  one  epoch  are  of  limited  use,  but  they  become  of  enor- 
mous value  when  combined  with  several  epochs  and  changes  in  the  positions  are 
observed.  From  studying  stellar  positional  changes  several  important  astronomi- 
cal parameters  can  be  determined;  probably  the  most  important  of  these  are  stellar 
parallaxes  and  proper  motions.  Trigonometric  parallaxes  are  the  only  method  for 
determining  stellar  distances  based  only  on  geometry.  Once  the  distances  of  stars 
are  known,  then  a  host  of  other  quantities  can  be  found,  for  example,  absolute 
magnitudes  and  in  double  star  systems,  actual  orbit  semi-major  axes  and  the  sum 
of  the  masses  of  the  stars  can  be  determined.  Proper  motions  in  combination  with 
radial  velocities  allow  one  to  study  the  motion  of  stars  in  clusters,  to  identify  stellar 
associations,  to  determine  the  orbital  motion  of  double  and  multiple  star  systems 
and  to  study  the  motions  within  the  Galaxy  which  leads  to  an  understanding  of 
the  evolution,  structure  and  rotational  curve  of  the  Galaxy. 

The  considerable  distances  to  stars  reduce  their  proper  motion  to  very  small 
amounts;  for  most  stars  these  are  considered  negligible  within  the  current  proce- 


3 
dures  of  measuring  precision.  It  is  almost  always  an  angle  too  small  to  measure 
with  precision  in  a  single  year.  The  majority  of  stellar  proper  motions  can  only 
be  determined  after  a  cumulative  effect  (as  a  rule)  over  many  decades  produces  a 
measurable  change  in  the  stars'  position.  Most  of  the  detectable  proper  motions 
are  for  nearby  stars.  Though  there  are  several  hundred  stars  with  proper  motions 
greater  than  l."0  per  year  (the  largest  motion  being  that  of  Barnard's  Star,  10. "25 
per  year),  the  mean  proper  motion  of  all  naked-eye  stars  is  less  than  0."1  (Abell 
et  al..  1987).  Therefore,  in  order  to  detect  accurately  proper  motions  we  must 
have  extremely  accurate  and  precise  positions  for  the  stars  with  preferably  decades 
having  elapsed  between  observations.  Star  catalogs,  in  particular  old  star  catalogs, 
are  extremely  useful  for  the  determination  of  proper  motions  for  a  large  number  of 
stars.  For  example,  stellar  proper  motions  can  be  detected  by  comparing  catalogs 
containing  the  same  stars  with  epochs  several  decades  apart  or  by  combining  the 
positions  found  in  a  single  catalog  with  new  observations  of  the  same  stars.  Thus, 
one  way  to  improve  upon  the  current  proper  motion  estimates  is  to  improve  upon 
the  accuracy  of  the  star  positions  found  in  the  old  star  catalogs  by  using  improved 
reduction  techniques  and  better  reference  material  than  was  originally  available. 
This  is  one  of  the  justifications  for  the  work  reported  on  in  this  paper. 

For  this  project  we  have  re-reduced  (re-determined)  the  star  positions  found 
in  the  Southern  Polar  Zone  of  the  Yale  Photographic  Star  Catalogue  (-70°  to 
-90°  declination,  Lii  1971).    Because  of  the  scarcity  of  southern  hemisphere 
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observatories  in  the  first  half  of  this  century,  there  are  few  accurate  old  (by 
this  we  mean  50  years  or  more)  star  catalogs  containing  southern  declination  star 
positions.  Thus,  this  region  of  the  Yale  Zone  Catalogues  is  particularly  interesting 
for  it  contains  some  of  the  oldest  accurate  observations  (epoch  1955-1956)  taken 
with  photographic  plates  of  the  southern  polar  zone.  In  addition,  the  photographic 
plates  that  were  used  to  create  this  catalog  were  photographed  so  that  they  heavily 
overlapped  (50%  in  right  ascension  and  50%  in  declination);  thus,  many  of  the 
same  stars  are  on  several  plates.  This  heavily  overlapped  pattern  is  particularly 
well  suited  to  our  reduction  technique  which  has  improved  the  accuracy  of  the 
stellar  positions  and  will  lead  to  more  accurate  and  precise  proper  motions  when 
calculated  from  the  positions  found  in  the  improved  catalog. 

The  Southern  Polar  Zone  of  the  Yale  Photographic  Star  Catalogue  is  one  of 
many  zones  which  make  up  the  Yale  Zone  Catalogues.  In  1913  the  director  of  the 
Yale  Observatory,  Frank  Schlesinger  ("father  of  modem  astrometry")  initiated  a 
project  to  photograph  the  entire  sky.  He  planned  to  re-observe  by  photography  the 
positions  of  the  Astronomische  Gesellschaft  Catalogue  (AGK  series)  stars  and  then 
by  comparing  the  older  AGK  positions  with  the  new  Yale  determinations  to  detect 
the  proper  motions  of  stars  to  the  ninth  magnitude.  Ida  Barney  collaborated  with 
Schlesinger  on  this  project  and  after  Schlesinger' s  death  in  1942,  continued  on  her 
own  until  1959.  Between  1925  and  1959  they  (with  very  few  other  colleagues) 
published  some  20  volumes  of  accurate  positions  and  attainable  proper  motions 
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for  nearly  146.000  stars  in  the  zone  from  -30°  to  +30°,  +50°  to  +60°  and  +85° 
to  +90°  (Hoffleit,   1962). 

Dirk  Brouwer.  director  of  the  Yale  Observatory  from  1941  to  until  his  death 
in  January  1966.  continued  Schlesinger's  project  and  extended  it  to  more  southern 
declinations:  -30°  to  -90°.  Early  in  1960,  Brouwer  obtained  a  2-year  contract 
from  the  U.S.  Army  Map  Service  to  re-observe  all  the  stars  in  the  Cape  Zone 
Catalogues  for  the  zones  -30°  to  -50°  and  -60°  to  -90°,  involving  some  70,000 
stars,  in  order  to  improve  the  positions  and  particularly  the  proper  motions  in 
those  zones.  Dorrit  Hoffleit  was  appointed  to  supervise  the  project,  starting  at 
-30°  and  working  southward  (Hoffleit,  1983). 

The  southern  catalogs,  like  the  northern  ones,  were  divided  into  zones  (ex- 
tended regions  of  the  sky  bounded  by  parallels  of  declination).  Except  for  the 
polar  caps,  each  zone  was  covered  by  photographic  plates  centered  at  the  same 
declination,  with  the  right  ascensions  of  the  centers  spaced  such  that  the  vertical 
center  line  of  a  plate  almost  coincided  with  the  vertical  edges  of  one  of  its  neigh- 
bors (i.e.  overlap  50%  in  right  ascension).  This  produced  a  double  coverage  of 
the  sky.  The  overlap  in  declination  was  at  most  only  a  few  degrees  and  often 
none.  The  Southern  Polar  Zone  has  a  much  stronger  overlap  pattern  than  the  other 
zones.  This  region  consists  of  4  belts  of  zones  with  a  declination  overlap  of  50% 
between  the  zones  and  an  overlap  of  50%  in  right  ascension  within  the  zones.  The 
Southern  Polar  zone  was  photographed  according  the  scheme  given  in  Table  1 . 
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Table 

1:  Yale  Zones 

Number  of  PI 

ates 

Center  Declination 

Sepai 

■ation  in  R.A. 

24 

-75 

ih  Qom 

20 

-80 

1  12 

16 

-85 

1  30 

4 

-90 

6  00 

The  strong  overlap  pattern  is  shown  the  polar  plots  in  Figure  1.  The  first 
plot  shows  the  different  declination  zones  the  center  representing  the  Southern 
Celestial  Pole  (-90°)  and  the  other  zones  are  represented  by  concentric  circles, 
moving  outward  from  the  pole:  -85°,  -80°,  -75°,  -70°  and  -65°.  The  first  plot 
is  also  only  3  plates,  the  next  4  plots  are  graphs  for  the  plates  in  each  zone  and 
the  last  plot  shows  all  the  plates  for  the  region.  The  number  of  plates  contained 
is  each  zone  given  in  Table  1. 

The  Yale  Zone  Catalogues  have  been  of  enormous  importance  to  positional 
astronomy  and  in  particular  the  determination  of  stellar  proper  motions.  Together 
they  contain  hundreds  of  thousands  of  stars.  Schlesinger  had  planned  to  eventually 
cover  the  whole  sky;  however,  this  aim  was  not  quite  achieved,  the  zones  -1-60° 
to  4-85°,  -1-30°  to  -f-50°  and  -50°  to  -60°  are  still  lacking.  The  region  we  are 
concerned  with  is  the  -70°  to  -90°  zone.  This  was  the  last  Yale  zone  to  be 
photographed  (1955-1956)  and  has  the  largest  star  density  and  magnitude  range 
of  all  the  Yale  zones. 
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Original  Reduction 

The  plate  material  for  the  Yale  Southern  Polar  Zone  Catalogue  consists  of 
64  photographic  plates  with  a  sky  coverage  of  1  l°x  1 1°  each.  These  plates  were 
exposed  between  September  15,  1955  to  August  7,  1956  at  Sydney,  Australia, 
measured  at  the  Yale  Observatory  between  1963  and  1968  and  the  reduction  was 
completed  in  1971.  The  reduction  of  the  plates  (i.e.  the  determination  of  the 
equatorial  positions  of  the  images  on  the  plates)  was  completed  in  1971  (Hoffleit, 
1971).  Figure  2  (top)  shows  the  stars  for  this  catalog  on  a  polar  plot.  The  bottom 
graph  in  Figure  2  shows  the  references  used  to  re-reduce  the  catalog. 

In  order  to  explain  the  reasoning  behind  our  desire  to  reduce  this  catalog, 
we  will  briefly  outline  the  procedure  for  reducing  plates.  The  details  of  this 
are  given  in  Chapter  4.  The  reduction  of  photographic  plates  requires  the  use 
of  a  set  of  reference  stars,  whose  stellar  images  are  on  the  plate  and  whose 
equatorial  coordinates  are  known  by  previous  investigation.  In  the  traditional 
method  of  reducing  plates,  estimates  of  certain  plate  parameters,  which  are 
numerically  different  and  characteristic  for  each  plate,  are  determined  from  least 
squares  reductions  of  the  reference  stars'  equatorial  coordinates  with  respect  to 
the  rectangular  coordinate  on  the  plate.  These  plate  parameters  are  then  used 
to  convert  field  stars  images  (i.e.,  nonreference  stars)  to  equatorial  coordinates. 
This  method  of  reducing  each  plate  separately  has  become  known  as  the  "single 


plate  method." 

The  original  reduction  of  the  Southern  Polar  Zone  was  performed  using  a 
single  plate  method  by  Phillip  Lu  and  Dorrit  Hoffleit.  The  reference  catalog  used 
was  the  Second  Cape  Catalog  for  1950.  The  Yale  plates  were  taken  between  1955 
and  1956,  so  the  reference  star  positions  taken  from  the  Cape  were  corrected  to 
the  epochs  of  the  Yale  plates  by  applying  proper  motions  for  intervals  of  three  to 
nine  years.  On  the  average  there  are  150  to  200  reference  stars  per  plate. 

After  the  reduction,  the  resulting  positions  of  stars  that  had  been  recorded  on 
more  than  one  plate  were  averaged  together.  The  resulting  catalog  consists  of 
18,702  stars  with  their  positions  and  proper  motions  on  the  FK3  system  (Third 
Fundamental  Katalog,  Astronomisches  Rechen-Institut,  1957).  "The  tabulated 
probable  errors  of  the  positions  represent  simply  the  internal  consistency  between 
different  Yale  plates  used  for  each  star.  They  imply  that  the  average  probable 
error  of  each  coordinate  determined  from  a  single  plate  amounts  to  ±0."33.  This 
is  consistent  with  the  average  formal  probable  error  of  a  plate  solution,  namely 
±0."35"  (Lu,  1971,  p.  29).  Note  that  the  (usually  quoted)  standard  errors  are 
about  1.5  times  this  amount  or  ±0."53. 

Why  Do  a  Re-Reduction? 

A  major  drawback  of  the  single  plate  method  is  that  it  produces  a  different 
position  estimate  (right  ascension,  declination)  for  each  star  on  each  plate,  even  if 
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that  star  also  appears  on  more  than  one  plate.  Clearly  the  star  can  occupy  only  one 
position  at  one  time,  but  this  cannot  be  obtained  without  mathematically  enforcing 
this  constraint  on  a  multiple  plate  process.  Instead,  the  separately  computed  values 
are  averaged  to  give  the  "best"  (in  the  mathematical  statistics  sense)  estimate  for 
the  stellar  position. 

As  just  mentioned,  the  most  important  constraint,  which  was  never  enforced, 
is  the  simple  fact  that  any  one  star,  at  a  given  instant  of  time,  occupies  only  one 
position.  In  the  case  of  the  Yale  Southern  Polar  Zone,  the  plates  were  so  strongly 
overlapped  that  most  stars  occurred  on  2  plates,  some  even  on  20  plates.  Thus  after 
performing  the  single  plate  reduction  there  were  as  many  as  20  different  position 
estimates  for  some  of  these  stars.  Taking  advantage  of  the  rich  overlap  we  have 
re-reduced  the  plates  using  a  more  powerful  reduction  technique,  the  overlapping 
plate  technique  (Eichhom,  1960;  for  review  Eichhom,  1985).  In  this  method  the 
critical  constraint  that  a  star  can  only  have  one  position  at  one  time  is  enforced. 
Then  all  the  reduction  parameters  on  all  plates  are  obtained  simultaneously  by 
solving  one  large  system  of  equations  (i.e.,  by  a  block  adjustment).  Whereas  the 
single  plate  is  limited  to  individual  least  square  reductions  of  just  the  reference 
stars  on  each  plate,  the  overlapping  plate  method  reduces  all  the  plates  at  the 
same  time  using  all  the  stars  (not  just  the  reference  stars).  More  accurate  (because 
the  systematic  errors  on  individual  plates  are  diminished  as  a  consequence  of  the 
forced  plate-to-plate  agreement  which  must  now  occur)  and  thus  also  more  precise 
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positions  than  those  obtained  from  the  single  plate  reduction  are  generated  by  the 
overlap  method.  This  method  also  has  many  capabilities  that  are  not  possible  with 
a  single  plate  reduction,  such  as  reducing  systematic  errors  between  neighboring 
plates  and  the  estimation  of  parameters  common  to  certain  data  sets  at  varying 
levels. 

Though  the  theoretical  formation  of  the  overlap  method  was  known  when  the 
original  reduction  was  done,  it  was  not  used.  Originally  it  had  been  planned  to 
use  an  overlap  reduction,  but  as  Hoffleit  explained  in  the  preface  of  this  catalog, 
"the  planned  more  elaborate  reductions  may  be  a  long  time  forthcoming  because 
of  a  present  shortage  of  both  funds  and  personnel"  (Hoffleit,  1971,  p.  3).  In  fact, 
when  the  computations  for  the  Yale  Catalogue  were  carried  out,  it  would  have 
been  very  difficult  to  take  advantage  of  the  huge  number  of  existing  constraints, 
because  enforcing  them  requires  a  computing  effort  that,  while  almost  modest  by 
contemporary  standards,  was  totally  prohibitive  at  the  time.  Many  of  the  existing 
geometric  and  physical  constraints  could  not  be  enforced  and  were  therefore 
wasted.  The  solution  using  all  the  available  information  has  become  possible 
only  through  the  emergence  of  fast  computers  with  large  memory  capacities. 

In  addition,  we  were  able  to  use  a  more  precise  and  accurate  reference  star 
catalog,  the  International  Reference  Stars  (IRS,  Corbin,  1991)  which  was  made 
available  to  us  by  the  US  Naval  Observatory  (see  Figure  2,  bottom).  Thus  by 
adjusting  the  originally  measured  position  of  the  stars'  images  on  the  plates  on 
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the  basis  of  better  reference  star  position  estimates  (which  were  not  available 
when  the  original  reductions  were  carried  out)  and  using  a  superior  adjustment 
algorithm,  the  early  positions  of  these  stars  have  been  improved  by  roughly  20%. 

The  improvement  in  the  positions  found  in  this  old  star  catalog  will  improve 
the  proper  motion  estimates  calculated  with  them.  For  example  the  ARCS  (Astro- 
graphic  Catalogue  Reference  Stars,  Corbin  and  Urban,  1991)  and  PPM  (Positions 
and  Proper  Motion  Catalog,  Bastian  et  al.,  1993)  catalogs  (both  key  reference 
catalogs)  used  the  Yale  Southern  Polar  Zone  Catalogue  for  their  determination  of 
proper  motions.  Thus  by  improving  the  positions  in  the  Yale  Catalogue  we  can 
improve  other  catalogs  that  use  the  positions  from  the  Yale  Catalogue.  This  is 
especially  important  in  the  Southern  Celestial  Hemisphere,  where  there  are  few 
old  reference  catalogs.  The  wholesale  improvement  of  the  precision  and  accuracy 
of  proper  motion  estimates  is  one  of  the  keys  to  improving  the  precision  and 
accuracy  of  a  large  variety  of  astronomical  data.  This  is  one  of  the  reasons  why 
the  measurements  in  the  Yale  Catalogue,  which  is  based  on  plates  taken  several 
decades  ago,  are  so  important. 


12 


Plates  1,  3  and  25 


■75°  Zone 


-80°  Zone 


-85°  Zone 


-90°  Zone  -70°  to  -90°  Zone 

Figure  1 :  The  plates  for  the  Southern  Polar  Zone  of  the  Yale  Catalogue 


13 


Star  Images  for  the  Yale  Southern  Zone 


IRS  Reference  Stars 

Figure  2:  The  Stars  of  the  Southern  Polar  Yale 
Zone  of  the  Yale  Photographic  Star  Catalogue 


CHAPTER  2 
PHOTOGRAPHIC  ASTROMETRY 

One  method  of  creating  a  star  catalog  (a  list  of  right  ascension  and  declination 
of  a  group  of  stars)  is  to  photograph  a  region  of  the  sky  with  photographic  plates. 
Once  the  plates  are  developed,  the  positions  of  the  star  images  on  the  plate  in 
a  certain  rectangular  coordinate  system  (x,y)  are  determined  with  a  measuring 
machine.  The  positions  of  the  measurements  on  the  photographic  plate  can  be 
derived  only  with  the  use  of  a  set  of  known  positions  of  stars  whose  images 
are  among  those  recorded  on  the  plate.  Their  equatorial  coordinates  are  usually 
found  from  existing  star  catalogs.  We  call  such  stars  reference  stars  and  all  other 
stars  for  which  we  do  not  have  equatorial  coordinates,  field  stars.  The  goal  of 
catalog  compilation  is  to  determine  accurately  the  spherical  coordinates  (for  our 
purposes,  equatorial  coordinates)  of  all  the  stars  imaged  on  the  plate  (field  stars 
as  well  as  reference  stars). 

In  the  ideal  case,  one  assumes  that  the  optical  system  (telescope  +  plate 
holder)  is  equivalent  to  a  pinhole  camera,  in  which  is  governed  by  the  laws  of 
the  gnomonic  projection.  A  gnomonic  projection  is  the  projection  of  a  spherical 
surface  onto  a  plane  through  a  point.  For  descriptive  purposes  we  will  assume  the 
existence  of  the  fictitious  "celestial  sphere"  and  we  will  further  assume  that  the 
center  of  the  celestial  sphere  lies  at  the  optical  center  of  the  objective  system. 
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In  this  perfect  case  the  portion  of  the  "celestial  sphere"  where  the  telescope 
is  pointing  is  projected  onto  a  plane  through  the  focal  point  of  the  telescope. 
The  resulting  rectangular  coordinates  projected  onto  the  plane  are  known  as 
the  standard  coordinates  i^j]).  The  usefulness  of  photography  for  astrometric 
positional  work  is  connected  with  the  fact  that  there  exists  a  simple,  rigorous  and 
parameter-free  geometrical  relationship  between  positions  of  objects  on  the  sky 
and  their  standard  coordinates  projected  onto  the  plate. 

Standard  Coordinates  and  their  Relationship  to  the  Equatorial  Coordinates 


tangent  plane 


celestial  sphere 


objective  =  center  of  celestial  sphere 
photographic  plate 


Figure  3;  Gnomonic  Projection 
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Imagine  a  plane  which  is  tangential  to  the  celestial  sphere  at  the  point  Qq,  ^o- 
In  the  plane  define  a  left  handed  coordinate  system  (<f,7/)  (as  seen  from  the  center 
of  the  sphere),  such  that  its  origin  is  at  the  tangential  point.  The  positive  ?;  axis 
is  tangential  to  the  hour  circle  through  Oq  and  points  northward.  The  (  axis  is 
tangential  to  the  parallel  of  declination  at  do  and  points  east  (in  the  direction  of 
increasing  right  ascension).  This  plane  coordinate  system  is  thus,  by  definition, 
parallel  to  the  focal  plane  of  the  telescope  and  therefore  parallel  to  the  plane  of 
the  photographic  plate  being  used  to  record  the  images.  The  tangential  point  can 
also  be  describesd  as  the  point  in  which  a  line  normal  to  the  plate  coincides  with 
the  optical  axis  of  the  telescope  and  penetrates  the  plate  close  to  its  geometrical 
center.  Thus  star  positions  on  the  celestial  sphere  are  projected  in  straight  lines 
through  the  center  of  the  sphere  (which  is  the  also  the  focal  point  of  the  telescope) 
and  onto  the  tangent  plate.  This  geometrical  situation  is  illustrated  in  figure  3. 

If  we  further  assume  that  the  radius  of  the  sphere  is  equal  to  1 ,  the  well  known 
relationships  between  a  point  on  the  celestial  sphere  (a,  6)  and  the  gnomonic 
projection  (i^,  rj)  of  this  point  onto  a  plane  which  is  tangent  to  the  unit  sphere 
at  iQo,6())  is  given  by 


cos  ^ sin  (q  —  ao) 

sin  6  sin  Sg  +  cos  8  cos  Sg  cos  (a  - 
sin  6  cos  So  —  cos  6  sin  Sg  cos  (o  - 

-Qo) 

-Olo) 

(1) 

sin  6  sin  Sg  +  COS  6  COS  6o  COS  {a  —  Og). 
These  equations  are  rigorous  and  free  of  estimated  parameters,  and  can  be  inverted 

so  that  Q  and  S  can  be  found  from  ^  and  //.   Note  that  ^  and  r]  are  in  units  of 
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the  telescope's  focal  length. 

However,  in  practice  no  real  optical  device  produces  images  in  the  same  way 
as  a  pinhole  camera.  One  needs  to  take  into  account  such  effects  as  coma,  radial 
distortion,  rotation  and  tilt  of  the  plate  when  photographed  or  measured.  The 
measured  coordinates  (x,y)  are  actually  approximations  to  the  standard  coordinates 

(s^^)• 

Determining  the  Equatorial  Coordinates  from  the  Measured  Coordinates 

The  goal  of  a  plate  reduction  is  the  derivation  of  the  spherical  equatorial 
coordinates  from  the  objects'  measured  positions  on  the  plate.  Usually  this  process 
is  achieved  in  two  steps.  First  the  measured  coordinates  are  transformed  to  the 
standard  coordinates,  then  these  are  converted  to  equatorial  coordinates.  The 
second  step  has  already  been  dealt  with  (i.e.,  equation  1),  so  only  the  first  remains. 

The  establishment  of  the  relationship  between  the  standard  coordinates  {^,t]) 
and  the  measured  plate  coordinates  (x,y)  of  the  star's  images  is  the  fundamental 
problem  of  determining  the  equatorial  coordinates.  In  the  ideal  case  (assuming  a 
pinhole  camera),  even  though  the  precise  position  of  the  image  of  the  tangential 
point  (qq,  6o)  on  the  plate  is  not  known,  it  is  reasonable  to  assume  that  it  is  near 
the  geometric  center  of  the  plate.  Plates  are  usually  exposed  in  such  a  manner 
that  it  is  also  reasonable  to  assume  that  the  edges  of  the  plate  (i.e.  the  x  and  y 
axes)  are  closely  parallel  to  the  ^  and  r;  axes,  respectively.   If  the  origin  of  the 
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measured  coordinates  is  chosen  to  lie  at  the  geometric  center  of  the  plate  and 
the  axes  are  properly  oriented,  to  first  approximation  the  standard  coordinates  are 
related  to  the  measured  coordinates  by  a  factor  of  the  telescope's  focal  length,  s: 

(2) 

The  unit  of  length  on  the  plate  in  the  ^-r]  system  is  such  that  the  standard 
coordinates  are  equal  to  1  at  a  distance  of  1  focal  length  from  the  tangential 
point.  Each  departure  from  the  ideal  imaging  process,  a  gnomonic  projection, 
results  in  deviations  from  the  perfect  relationship  in  the  equation 

f"-^U.f^°^^       ''^^^](^)  (3) 

\y  —  a  J  y— sm0     cos  0  J  \r]  J 

where  c,  d  and  <p  are  all  small  quantities  of  the  first  order.  The  end  result  is  a 
mapping  between  the  standard  coordinates  and  the  measured  coordinates.  Tradi- 
tionally, this  mapping  is  called  a  plate  model.  The  plate  parameters  have  tradition- 
ally been  called  plate  constants.  We  prefer  the  terminology  "plate  parameters" 
which  indicates  that  these  are  adjustment  parameters  to  be  estimated  by  some 
statistical  process. 

In  the  traditional  method  of  reducing  plates  one  assumes  there  exists  for  each 
plate  the  following  relationships 

(4) 
g(xi,y,;mi,  c,;a,,^,;ai.a2,  ...,a„)  =  0 
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between  the  measured  coordinates  x,,  y,  of  the  image  of  the  /'th  star,  its  magnitude 
An,  and  color  index  q,  and  its  coordinates  o„  6 a,  on  the  one  hand  and  a  set  of  plate 
parameters  {a^ }  on  the  other  hand.  If  all  the  plates  had  been  taken  with  the  same 
telescope,  the  plate  parameters  will  belong  to  groups  such  that  some  elements  vary 
from  plate  to  plate,  others  vary  from  plate  to  plate  but  are  constrained,  and  finally 
some  are  the  same  for  all  plates  obtained  on  the  same  telescope.  (Eichhorn,  1974). 

After  the  functional  form  of  the  plate  model  is  established,  the  plate  constants 
(or  estimates  thereof)  are  determined  from  a  statistical  adjustment  (usually  a  least 
squares  adjustment),  by  fitting  the  measured  coordinates  to  the  corresponding 
spherical  coordinates  of  the  reference  stars  (a,  6).  In  this  case,  equations  (4)  then 
become  the  equations  of  condition  for  a  least  squares  adjustment  in  which  the  x,, 
y,,  Q,  and  (5,  are  regarded  as  observations  and  the  a/,  a2,  ...,  a„  as  the  unknown 
parameters.  The  m,  and  c,  are  assumed  to  be  known  for  every  star.  Of  course, 
there  must  be  enough  reference  stars  to  make  the  system  over-determined  (i.e., 
the  number  of  reference  stars  must  be  greater  than  the  number  of  parameters).  A 
good  rule  of  thumb  is  that  there  should  be  at  least  3  times  as  many  reference  stars 
as  plate  parameters.  Once  the  plate  parameters  are  determined  the  best  estimates 
of  the  standard  coordinates  for  all  the  stars  (field  stars  as  well  are  reference  stars) 
can  be  obtained.  Then,  after  routinely  inverting  the  gnomonic  projection  (equation 
1)  the  equatorial  coordinates  for  all  the  stars  are  computed. 
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The  Relationship  Between  the  Standard  and  Measured  Coordinates 

Even  under  the  most  favorable  conditions  the  assumption  of  a  perfect 
gnomonic  projection  can  only  approximate  the  actual  circumstances.  It  is  clear 
that  neither  telescope  nor  measuring  machine  can  be  constructed  with  perfect  ac- 
curacy. First  of  all,  there  is  no  optical  system  in  practical  use  which  unites  in  a 
single  point  all  light  rays  that  enter  the  front  component  of  the  objective  from  the 
same  direction.  In  reality  there  are  many  deviations  from  the  ideal  situation  which 
require  small  corrections  to  be  added  to  equation  (2).  These  departures  from  the 
perfect  gnomonic  projection  are  "instrumental."  Others,  which  rotate  and  distort 
the  coordinate  system  with  respect  to  which  the  position  estimates  of  the  reference 
stars  are  catalogued,  may  be  termed  astronomical  effects  and  include  refraction, 
aberration,  precession  and  nutation.  We  will  deal  separately  with  these  two  types 
of  "corrections."  The  "instrumental"  corrections  will  be  characterized  by  the  plate 
model  and  the  "spherical"  corrections  will  be  discussed  in  the  next  section. 

Below  is  a  list  of  some  of  the  more  common  departures  from  the  ideal 
gnomonic  projection;  this  list  is  by  no  means  exhaustive: 

1.  Shift  in  the  origin  of  the  x-y  coordinate  system  with  respect  to  the  ^-rj 
coordinate  system  (c  and  d  in  equation  3 

2.  Rotation  of  the  ^-j]  axis  with  respect  to  the  x-y  axis. 
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3.  The  ^-j]  scale  may  be  different  from  the  x-y  scale,  due  to  a  poorly  known 
focal  length. 

4.  Tilt  of  the  x-y  focal  plane  with  respect  to  the  ^-r]  focal  plane,  this  will  cause 
an  incorrect  tangential  point. 

5.  Different  scales  on  each  axis,  due  to  different  measuring  screws  or  otherwise 
introduced  by  the  measuring  device. 

6.  Lens  aberrations  (distortion,  spherical  aberration,  coma,  astigmatism,  chro- 
matic aberration)  affect  the  position  of  the  stellar  centroids  either  by  enlarging 
them  (thus  making  it  more  difficult  to  determine  the  center)  or  by  an  actual 
shift  in  positions,  both  of  which  vary  depending  on  the  location  on  the  plate. 
These  effects  plus  peculiarities  of  the  process  of  photographic  image  forma- 
tion render  the  positions  of  the  star's  images  on  the  plate  dependent  on  their 
magnitudes  and  spectra  expressed  in  terms  of  color  indices.  These  conditions 
lead  to  the  inclusion  of  the  m,  and  c,  terms  in  the  mapping  between  x-y  and 
^-T}  coordinate  systems. 

7.  Besides  lens  aberration  "smeared  out"  or  "blow  up  "  images  can  be  produced 
by  guiding  errors,  developing  errors  (e.g.  emulsion  shifts),  and  bad  seeing. 
These  can  and  do  affect  the  position  of  the  stellar  centroids. 

Each  departure  from  the  ideal  imaging  process  results  in  deviations  from  the 
relationship  x=s^\  y=57/.  Given  below  are  the  relationships  between  the  standard 
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coordinates  and  the  measured  coordinates  to  correct  for  many  of  the  commonly 
experienced  errors. 

1.  Translation  error:  If  the  photographic  plate  is  not  centered  correctly  on 
the  measuring  machine,  there  will  be  a  systematic  difference  between  the 
measured  coordinates  and  ideal  coordinates. 


^^- 


(c,d) 


Figure  4:  Origin  Shift 


X  —  s(  =  c 
y-  srj  =  d 


(5) 


23 
2.    Rotation:  If  the  photographic  plate  is  correctly  centered,  but  rotated  relative 

to  the  measuring  machine  axes  by  an  angle  (p  {(p  >  0  for  counterclockwise). 

X  =  s(^  cos  6  +  r]  sin  (p) 

(6) 

y  =  .s(;;cos  c!)  —  .^sin  0). 


T1 


U 


.^ 


Figure  5:  Rotation  of  x-y  axis  with  respect  to  the  ^-77  axis 


3.    Imprecise  focal  length:    An  incorrect  focal  length  can  cause  differences  in 
the  scales  between  the  ^,r]  system  and  x,y  system.  This  can  be  corrected  by 


making  small  changes  in  if  and  r/: 

X  —  Sif  =  a,f 

y  -  st;  =  er/. 


(7) 
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4.  Center  error  (tangent  point  error):  For  wide  fields  (  say  over  5°  x  5°)  the 

extraneous  methods  for  the  determination  of  qq  and  6o  yield  approximate 

values  of  these  equations  and  thus  they  must  be  corrected  differentially.   If 

the  actual  position  of  the  tangential  point  is  (oo+dao),  (^o+d<^o)  instead  of 

cvo,6o,  values  of  ^,i]  will  be  incorrect  if  calculated  only  from  qo,^o-  It  can  be 

shown  to  first  order  in  doo  and  d^o  (Eichhorn,  1971)  that  the  centering  error 

can  be  corrected  (to  first  order)  by 

X  —  s^  =  (cos  Sodoo)^^  +  dSo^T] 

(8) 
y  -  srj  =(cos  6odao)^T}  +  dSoT]^ 

provided  that  the  plate  was  oriented  during  measurement  so  that  the  x  and  ^ 
axis  as  well  as  the  y  and  //  axes  are  nearly  parallel  and  (^,  rj)  are  in  radians 

If  doo  and  d^o  turn  out  to  be  large,  it  is  recommended  that  the  calculation  be 
repeated  with  the  newly  calculated  (,,ri  now  referred  to  the  tangential  point  oo+dao 
and  ^o+d^Q.  This  may  be  repeated  until  further  repetitions  result  in  no  change  in 
the  assumed  qq,  6q  (Eichhorn,  1974). 

5.  Tilt:  If  the  photographic  plate  is  not  perpendicular  to  the  optical  axis  of 
the  telescope  but  titled  by  an  angle  u  relative  to  the  focal  plane,  the  errors 
introduced  will  be 

y  -  ST]  =  [p^T]  +  qr?-)  tana; 
where  p  and  q  are  constants  (Taff  1981).    One  can  see  the  terms  for  the 

corrections  for  zero  point  error  and  tilt  have  the  same  form.  This  is  because  if 
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we  have  a  situation  where  there  is  a  tangential  point  error  this  is  the  equivalent 
of  having  an  imaginary  plate  whose  true  tangent  point  coincides  with  our 
assumed  one,  but  is  inclined  to  the  actual  plate.  Since  these  two  effects  can 
be  modeled  in  a  similar  manner,  they  are  sometimes  combined  under  that 
same  name.  If  construction  and  proper  alignment  of  the  telescope  has  been 
carefully  executed  such  second  order  errors  may  be  neglected.  However,  in 
the  case  of  large  plates  (over  5°  x  5°),  it  is  important  that  these  higher  order 
terms  be  included  (Konig,  1962). 
6.  Radial  distortion:  If  the  image  process  produces  a  net  radial  distortion,  it  can 
be  modeled.  Radial  distortion  originates  at  the  intersection  point  of  the  optical 
axis  with  the  focal  plane  and  its  origin  will,  therefore,  also  be  very  close  to  the 
tangential  point.  This  aberration  is  radially  symmetric  and  usually  modeled 
in  the  form 

n 


'~n  (10) 


k=\ 
where  R  is  a  constant  (Eichhom  1974).   Unless  the  radial  distortion  is  very 

strong,  and  for  the  field  sizes  used  on  photographic  catalog  astrometry,  only 

the  first  terms  of  equations  (9)  are  necessary  for  appropriate  modeling  of 

radial  distortion.   Thus  with  n  =  1, 

(11) 
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7.    Nonperpendicularity  of  axes:  If  the  axes  of  the  measuring  machine  are  not 
perpendicular  to  each  other  (in  other  words  x  is  not  perpendicular  to  y),  let 


V  be  the  angle  between  the  r)  and  y  axes  (McNally,  1975) 

X  —  S(^=;/  tan  I/!' 

y  —  S77=(l  —  sec  t/')r/. 


(12) 


Figure  6:  Nonperpendicularity  of  axes 


8.    Coma:    The  effective  focal  length  of  the  instrument  may  depend  on  the 
apparent  magnitude  m  of  the  object  that  is  being  imaged.    This  is  coma; 


it  is  radial  and  accounted  for  by 

X  —  s^=tmx 

y  -  S7]=tmy. 


(13) 
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Experience  shows  that  the  coefficient  s  is  rather  sensitive  to  the  rate  of  change  of 
the  ambient  temperature  during  the  exposure.  Thus  it  can  not  be  regarded  as  a 
constant  for  any  one  objective,  although  the  values  obtained  in  practice  will  have 
a  tendency  to  cluster  around  a  certain  value  (see  Eichhom  &  Gatewood,  1967). 
9.    Decentering  distortion:  If  all  the  components  of  the  objective  are  not  prop- 
erly aligned,  the  resulting  imperfection  is  known  as  decentering  distortion. 
Brown  (1966)  found  that  the  appropriate  corrections  to  the  measured  x  and 
y  coordinates  are 


i+f^p,,+,{e+ri'y 


k=i 

n 


(14) 


k=\ 
again,  under  the  conditions  imposed  by  catalog  photographic  astrometry,  the 

terms  with  Pk+2    may  usually  all  be  assumed  to  equal  zero,  this  results  in  a 

second  order  effect 

(15) 
y  -  srj=2[P,^ri+P,ri']  +  P^i^+ri^)- 
One  can  see  by  comparison,  the  first  brackets  are  exactly  the  same  structure  as 

those  that  model  the  correction  to  the  tangential  point  error  (number  4)  and  tilt 

(number  5).    If  the  model  includes  the  provisions  for  correcting  the  tangential 

point  (or  tilt)  anyway,  the  above  decentering  terms  in  the  brackets  can  be  may 

"lumped"  with  them.   Then  by  just  adding  the  Pi(^^  +  r/^);  P2(^^  +  7;-)terms 

one  can  account  for  the  corrections  due  to  decentering  distortion. 


28 
As  just  indicated  the  sources  of  errors  are  not  independent  of  each  other 
and  there  is  much  crossover  in  effects  but  they  can  be  modeled  in  a  first  order 
(sometimes  second  order)  approximation  that  will  more  than  suffice  considering 
their  magnitude  and  the  magnitude  of  the  measuring  errors.  It  is  worth  noting 
that  the  coefficients  which  describe  the  distortions  are  very  insensitive  to  changes 
in  the  ambient  temperature  and  may  thus  be  regarded  as  more  or  less  constant 
for  any  one  objective.  This  statement  is  not  true  for  coma,  which  is,  in  effect,  a 
dependence  of  the  effective  focal  length  (scale)  on  magnitude. 

Some  Commonly  Used  Models 

The  difficulty  in  establishing  the  relationship  between  the  standard  and  mea- 
sured coordinate  system  is  in  determining  which  terms  to  include  in  the  model. 
Eichhom  (1985)  writes  the  relationship  (2)  in  the  form 

where  S,  the  model  matrix,  is  a  function  of  (,  t],  m  and  c,  and  a  is  the  vector 
of  plate  parameters.  The  form  of  the  model  matrix  is  determined  by  probing 
possible  deviations  from  the  gnomonic  projection,  first  adding  appropriate  terms 
to  a  very  basic  model,  and  then  investigating  the  adjustment  residuals  produced  by 
additional  terms.  When  significant  terms  are  found  they  are  added  to  the  model. 
This  trial  and  error  process  continues  until  one  is  satisfied  that  one  has  found  the 
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"best"  possible  model  which  approximates  the  "true"  system.  We  will  now  look 
at  some  of  the  models  which  have  been  employed  in  the  past. 

4-Constant  Model 


The  first  model  we  will  look  at  is  the  simplest  and  therefore  the  most  basic 
model;  it  corrects  for  errors  in  the  measurements.  We  will  assume  the  tangential 
(ao  ^o)  to  be  known  or  determinable  by  extraneous  methods  and  thus  not  subject 
to  the  adjustment.  The  scale  and  orientation  of  the  standard  coordinates  are  fixed 
by  their  definition.  However,  the  system  of  the  measured  coordinates  is  dependent 
on  the  measuring  machine  and  the  accidental  orientation  of  the  photographic  plate 
on  it.  Consequently,  the  two  systems  may  not  agree  because  the  origins  of  each 
system  might  be  different  and  the  axes  of  one  system  may  be  rotated  with  respect 
to  the  ones  of  the  other.  This  situation  is  demonstrated  in  figure  7. 


y 

V 


1 
1 

1 

^_._._ 

T(C.D) 

Figure  7:  4-Constant  Model:  origin  shift  and  rotation  of  axes 


(19) 
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X  =  s{^  cos  d  +  t]  sin  o)  +  c 

(17) 
y  =  .1(7;  cos  6  —  ^s'md))  +  d. 

where  o  is  the  angle  of  rotation,  s  the  focal  length  and  c  and  d  the  difference  in 

the  origins.  Combining  these  effects  results  in  the  following  model.  If  we  let  a 

=  s  coso  and  b  =  s  sino,  then  the  equation  becomes, 

\  =  a^  -\-  br}  +  c 

(18) 

y  =  ar/  —  6(^  +  d. 

Frequently,  equation  (18)  is  written  in  the  form 

/a\ 

'^     n     1    o\    b 

,1       -e      0       lyl       C 

In  this  case  plate  model,  E,  is  the  2x4  matrix  and  the  vector  of  plate  constants 
a  is  (a,  b,  c,  and  d).  There  are  four  plate  constants,  so  this  model  is  called  the 
4-  constant  model. 

Usually,  the  plates  are  measured  in  such  a  way  that  the  x-axis  is  almost 

parallel  to  the  i^-axis  which  will  (if  x  and  y  are  strictly  rectangular)  render  the 

y-axis  almost  parallel  to  the  7;-axis.  This  will  make  b  very  small  compared  to 

a.    So,  when  4)"^  0, 

a  =  pcos  (f>        cos<^  Si  1, 

(20) 

thus,  a  ^  p 

and 

6  =  /9sin^       sin  0  si  0, 

(21) 

therefore,  6  si  0. 
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Thus  the  term  a  is  the  correction  to  the  assumed  focal  length.  If  the  x-y  and  ^-i] 
axes  are  tilted  with  respect  to  each  other,  this  term  may  have  a  small  effect  due 
to  this  rotation,  but  the  dominant  effect  is  the  correction  to  the  focal  length. 
The  inversion  of  the  above  equation  is 


1 


a^  +  6M^ 


(22) 


or  after  introducing  certain  quantities 


i\  A      B 

,,)       \-B     A 


iH'o 


(23) 


where  A,  B,  C,  D  make  up  the  matrix  of  plate  constants.  Thus  for  the  4-constant 
model  a  rigorous  equation  exists  for  ^  and  r/  to  be  solved  in  terms  of  the  plate 
constants  and  x  and  y. 

6-Constant  Model 

The  4-constant  model  does  not  correct  for  different  scales  in  the  x  and  y 
measurements.  Correcting  for  this  effect,  equation  (17)  becomes 


^  /  cos  (j)      sin  (?i>  \  /  A       /c 
'^y— sin0     cos<f)J\T]J       Id 


(24) 


or  after  multiplying  out 

X  —  s  cos  (l>^  +  s  sin  (f>Tj  +  c  —  a^ 

y  =  —s  sin  (f)^  +  s  cos  <f)r]  +  d  —  erj. 


(25) 
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After  substituting  the  following  values,  a  =  cos  6  -  ^,  e  =  cos  o  -  -,  6  =  sin  6, 

/"  =  -  sin  o,  c=-.  and  d  =  - 


sM    ,   /s^     '/     1     0     0     0 
'  0     0     0     ^      r;      1 


b 
c 
e 
f 


(26) 


1.  a  &.  e:  correction  to  assumed  focal  length. 

2.  b  &.f:  correction  for  rotation  of  x-y  system  with  respect  to  ^-r]  system.  This 
also  corrects  for  nonperpendicularity  of  the  x  and  y  axes.  These  should  be 
approximately  equal  and  opposite  in  sign. 

3.  c  &  d:  differences  in  origin  of  the  x-y  system  and  ^-q  system. 

For  the  6-constant  model  there  exists  an  equation  analogous  to  equation  (23) 
for  the  4-constant  model.  In  other  words  ^  and  r]  can  be  rigorously  solved  for 
from  the  plate  constants  and  x  and  y.  This  is  not  true  for  the  other  models  we 
will  be  using.  In  those  cases,  an  initial  approximation  is  used  and  the  solution 
then  iterated  until  a  set  tolerance  is  met. 

12-Constant  Model 

So  far  our  models  have  only  dealt  with  terms  linear  in  the  standard  coordinates 
(f  and  rj.  However,  there  are  many  "second  order"  effects  that  need  to  be 
considered,  especially  when  dealing  with  large  plates  like  the  Yale  plates  (11° 
X  11°).  A  standard  model  for  the  conversion  of  the  measured  coordinates  to  the 
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standard  coordinates  considers  corrections  for  I)  rotation,  2)  translational  error, 
3)  tangential  point  correction,  4)  radial  distortion,  5)  tilt  and  6)  coma  ;  it  is  the 
12-constant  model: 


+ 


0  0 


1  0  0 

0  ^  n 


0  C-  ^rj  m-mo        0        fm  ^{e  +  V') 

1  ^T]  J]-         0         m  -  nio  r?m  r]{^- +  rj"^) 
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f 

e 
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(27) 


1 .  a  &  e:  correction  to  assumed  focal  length  (different  in  x  and  y). 

2.  b  &  f:    correction  for  rotation  of  x-y  system  with  ^-tj  system.    This  also 
corrects  for  nonperpendicularity  of  the  x  and  y  axes 

3.  c  &  d:  differences  in  origin  of  the  x-y  system  and  ^-q  system 

4.  p  &  q:  Tangential  point  correction  and  tilt 

5.  /  &  j:   magnitude  "equations" 

6.  g:    coma 

7.  h:    radial  distortion 
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All  terms  of  Ha  are  at  least  one  order  of  magnitude  smaller  than  x,  y  or  ^,  q  and 
.V  is.  of  course,  an  approximation  for  the  focal  length.  This  clearly  shows  that  the 
rigorous  formulas  give  i^,  ?/  not  in  radians,  but  rather  in  units  of  the  focal  length; 
they  are  indeed  dimensionless  quantities. 

This  model  is  very  similar  to  one  we  chose  to  best  represent  the  spherical 
coordinate  system.  The  actual  model  and  method  of  determining  the  model  for 
the  Yale  plates  is  given  in  Chapter  6. 

Spherical  Correction 

We  stated  earlier  that  we  split  the  deviations  from  the  perfect  gnomonic 
projection  into  instrumental  and  astronomical  effects.  The  instrumental  effects 
are  corrected  by  the  plate  model.  The  astronomical  deviations  arise  because 
the  image  formed  on  the  photographic  plate  represents  the  refracted  apparent 
topocentric  positions  of  the  objects.  The  standard  coordinates  are  computed  from 
the  coordinates  listed  in  star  catalogs  (e.g.  mean  positions).  To  account  for  the 
difference  between  the  refracted  topocentric  apparent  positions  (x,  y)  and  the  mean 
positions  {^,  t}),  corrections  for  refraction  ,  precession,  and  for  extreme  accuracy, 
diurnal  aberration,  should  be  applied.  We  will  not  consider  diurnal  aberration, 
since  its  effect  is  smaller  than  our  measuring  error.  The  effects  of  precession  are 
accounted  for  by  correcting  the  reference  catalog  and  influence  the  results  only 
slightly  through  refraction. 
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Atmospheric  refraction  causes  light  from  a  celestial  body  to  be  bent  as  it  passes 
through  the  Earth's  atmosphere.  Assuming  a  simple  model  for  the  atmosphere  as 
a  plane  atmosphere  and  for  small  zenith  distances,  total  refraction  in  the  zenith 
distance  can  be  calculated  by 

z  —  zq  =  Az  =  /i  tan  z  (28) 

where  z  and  Zq  are  the  topocentric  apparent  (i.e.,  unaffected  by  refraction)  and 
observed  refracted  zenith  (i.e.,  effected  by  refraction)  distance  of  the  object, 
respectively.  3  is  the  constant  of  refraction  and  depends  weakly  on  the  temperature 
and  pressure  at  the  observing  station  (usually  this  value  is  about  1  minute  of  arc). 
Since  Az  >  0,  astronomical  refraction  raises  celestial  objects  toward  the  zenith. 
Astronomical  refraction  ideally  does  not  affect  the  azimuth. 

For  large  zenith  distances,  allowance  has  to  be  made  for  the  variation  of  3 
within  the  field  of  the  plate;  in  practice,  it  suffices  to  put  3  =  3o  +  i^'  tan'  z, 
where  3'  and  /3o  are  refraction  constants  and  can  be  found  in  refraction  tables 
(Konig,  1933).   Thus  equation  (28)  becomes 

Az  =  /3o  tan  z  +  /?'  tan^  z.  (29) 

We  have  neglected  higher  order  terms  for  zenith  distance  less  than  70°. 

We  will  split  the  effects  of  atmospheric  refraction  into  those  of  "absolute" 
and  "differential."  By  "absolute"  we  mean  the  amount  of  refraction  affecting  the 
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center  of  plate  or  more  precisely  the  tangential  point.  We  can  think  of  refraction 
as  shifting  all  the  positions  on  the  plate  by  this  constant  amount  plus  a  differential 
amount  which  varies  across  the  plate.  For  example  in  Figure  8  (top)  we  have 
plotted  the  effects  of  refraction  for  plate  1.  As  one  can  see  the  positions  are 
shifted  mostly  in  the  y  direction  (an  average  of  about  49").  In  Figure  8  (bottom) 
we  have  shown  the  effect  across  the  plate  by  subtracting  the  tangential  point 
refraction.  Note  that  the  scales  for  the  two  graphs  are  not  the  same;  in  the  bottom 
left  hand  comer  of  each  graph  we  have  plotted  the  scale:  50"  for  the  top  graph 
and  5"  for  the  bottom  one.  We  will  apply  a  correction  for  differential  refraction, 
and  since  our  plate  model  already  has  a  linear  term  in  it  we  will  allow  the  linear 
effect  of  absolute  refraction  to  be  taken  care  of  by  our  plate  model. 

It  has  been  a  frequent  practice  to  account  for  differential  refraction  by  adding 
higher  order  terms  to  the  model.  This  is  usually  achieved  by  adopting  a  full 
second  and  third  order  polynomial  expansion  in  the  coordinates.  There  are  two 
major  disadvantages  to  using  this  approach. 

Firstly,  if  the  number  of  parameters  is  large  an  unfavorable  error  propagation 
of  the  resulting  star  positions  in  the  plate  field  due  to  the  random  errors  in 
the  x-y  measurements  and  the  reference  star  positions  is  introduced  (Eichhom 
and  Williams,  1963);  in  particular  when  the  number  of  reference  stars  per  plate 
is  small.  Secondly,  all  effects  are  mixed  up  in  the  model  and  therefore  it  is 
impossible  to  separate  different  contributions  and  evaluate  the  physical  meaning 
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of  the  plate  parameters  in  detail. 

As  we  mentioned  above  we  will  pre-correct  the  x,y  measurements  for  the 
differential  atmospheric  refraction.  These  values  are  used  as  observations  in  the 
plate  reduction.  In  this  way,  the  number  of  parameters  in  the  model  can  be  kept 
to  a  minium  which  will  give  us  a  better  chance  to  disentangle  the  complicated 
origins  of  many  of  the  plate  parameters  (Zacharias  et  al.,  1992). 

The  usual  approach  in  plate  reduction  calculations  is  to  transform  the  equa- 
torial system  to  the  horizon  system  and  make  refraction  corrections  in  a  simple 
way  in  the  horizon  system.  In  this  manner  the  effect  of  refraction  is  added  to 
the  reference  positions.  However,  this  leads  to  the  situation  where  the  same  star 
imaged  on  different  plates  has  different  positions  due  to  refraction.  From  Chapter 
1  we  stated  that  the  power  of  the  reduction  technique  that  we  use,  the  overlap- 
ping plate  method  (Chapter  4),  is  that  it  enforces  the  constraint  that  a  single  star 
photographed  on  many  plates  must  have  only  one  position  in  the  sky.  Therefore 
correcting  the  stars'  equatorial  coordinates  for  refraction  leads  to  a  different  posi- 
tion depending  on  what  plate  it  was  imaged,  thereby  invalidating  our  constraint. 
We  have  chosen,  instead,  to  pre-correct  the  measured  star  images  so  that  they 
represent  the  position  on  a  plate  if  the  plate  had  been  observed  outside  of  the 
atmosphere.  We  have  used  the  equations  developed  by  Konig  (1933)  for  calcu- 
lating the  corrections  for  refraction,  Ax  and  Ay,  which  are  to  be  applied  to  the 
apparent  topocentric  coordinates,  x  and  y.  After  calculating  the  zenith  distance  of 


38 

the  plates  at  the  time  of  exposure  and  then  allowing  for  temperature  and  pressure 
we  can  determine  Jq  froni  the  refraction  Tables  7-13  in  Konig  (1933)  and  let  3' 
=  -0.082".  Table  2  and  Table  3  contain  all  the  necessary  information  to  determine 
refraction  with  the  aid  of  Konig' s  refraction  tables.  In  addition,  the  latitude  at  the 
Sydney  observatory  is  -33°  51'  41."1,  the  longitude  is  -10*^  4*"  49.'06,  and  we 
have  assumed  an  altitude  at  sea  level  and  a  barometer  reading  of  750  millibars. 
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Refraction  Correction:  Absolute 
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Refraction  Correction:   Relative 
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Figure  8:  The  effects  of  refraction  on  plate  1 
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The  tangential  point  is  taken  to  be  the  refracted  oq,  Sq.  The  rigorous  mathe- 
matical relations  between  the  refracted  positions  and  the  topocentric  positions  are 
somewhat  complicated.  For  practical  application,  only  development  into  series 
can  be  considered;  it  is  therefore  necessary  to  state  which  terms  must  be  retained 
and  which  can  be  neglected. 

In  order  to  obtain  the  simple  expressions  for  the  coefficients  of  the  series, 
Konig  introduces  the  following  auxiliary  quantities 

ki  =  tan  z  sin  \ 

ko  =  tan  z  cos  x 

(30) 

k3  =  l  +  kf 

A;4  =  1  +  ko 
where  x  is  the  (refracted)  parallactic  angle  at  the  tangent  point  (taken  to  be  the 
refracted  tangential  point)  and  ki    and  k2  have  a  simple  geometrical  meaning; 
they  are  the  tangential  coordinates  of  the  zenith  in  the  plane  of  the  plate.   The 
corrections  for  the  measured  values  are  (to  the  third  order): 
Ax=  [/?k3  +  2/?'(l  +  ki  +  ki)k2]x 

-  [/?ki(k2  +  UnSo)  +  2/?'(l  +  kj  +  k^)kik2]y 

+  2[/5  +  2/?'(l  +kj  +  ki)]kik2y  -  /?kik3x2  (31) 

-  2/9kfk2xy  -  /?kik4y2  +  f^kjx^  +  3/?kik2k3x2y 
+  .3{k3k4  +  2k^k^)xy^  +  /5kik2k4y^ 


Let 


and 


(33) 


(34) 
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Ay  =  [,.iki(k2  +  tan<5o)  +  2,J'(l  +  kj  +  k:])kik2]x 

+  [jk3  +  2i3'{l+ki  +  kl)kl]y 

+  [i  +  2J'(l+ki  +  k^)](k^-k^)y  (32) 

-  Jk2k3X-  -  2Jkikixy  -  jSk^Uy^  +  i^kikokax'^ 

+  /3(k3k4  +  2k2k^)x"y  +  3/:/kik2k4xy'  +  ^ik^vl 

Xr  =  -  l3kikix-  -  2/?kfk2xy  -  i3kiUy'  +  /3klx^ 

+3/^kik2k3x''y  +  /i(k3k4  +  2kik^)xy-  +  l^k^k-ik^y^ 

Yr  =  —  iSkyk^x"  —  'IlSkikixy  —  Jk2k4y"  +  /ikik2k3x"' 

+/?(k3k4  +  2k2k^)x2y  +  3/3kik2k4xy2  +  /5k|y3 
also  define 

A  =  /^k3  +  2/?'(l  +  ki  +  kyki 

B  =  -  [3ki{k2  +  tan)  +  2/?'(l  +  kj  +  k^kik.] 

C=2[,/?  +  2/:/'(l  +  kl  +  k^)]kik2y 

D  =  [^  +  2IS'{1  +kl  +  k',)]{kl-ki). 
Substituting  these  variables  into  equation  we  can  rewrite  equations  (31)  and  (32)  as 

Ax  =  Ax  +  By  +  C  +  Xr 

(36) 
Ay  =  -Bx  +  Ay  +  D  +  Yr. 

Written  in  this  form  the  refraction  terms  which  cause  a  change  in  the  scale  (A), 

rotation  (B),  linear  refraction  terms  (C  and  D),  residual  refraction  (XRand  Yr)  cab 

be  separated.  Generally,  for  small  fields  and  for  zenith  distance  less  than  70°  the 


(35) 
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terms  A  and  B  can  be  neglected.  However,  this  is  not  so  for  plates  taken  at  very 
high  declinations.  We  have  included  the  terms  for  change  of  scale  and  residual 
refraction,  but  not  rotation.  Rotation  was  not  included  because  the  rotation  caused 
by  refraction  cannot  be  separated  from  that  caused  by  the  position  of  the  plates 
in  the  measuring  machine. 
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Table  2:  Refraction 


PL 

RAO 

LMST 

MON 

1 

0.0 

0.108 

10 

2 

1.0 

1.042 

11 

3 

2.0 

2.058 

12 

4 

3.0 

2.75 

11 

5 

40 

4.342 

6 

5.0 

5.408 

7 

6.0 

6.158 

8 

7.0 

7.092 

9 

8.0 

8.092 

10 

9.0 

9.108 

11 

10.0 

10.063 

2 

12 

11.0 

11.258 

4 

13 

12.0 

11.958 

4 

14 

13.0 

13.083 

7 

15 

140 

14008 

4 

16 

15.0 

15.233 

6 

17 

16.0 

16.250 

6 

18 

17.0 

17.150 

6 

19 

18.0 

18.150 

6 

20 

19.0 

18.883 

6 

21 

20.0 

20.238 

9 

22 

21.0 

21.175 

10 

23 

22.0 

22.208 

10 

24 

23.0 

23.275 

10 

25 

0.0 

0.096 

11 

26 

1.2 

1.413 

10 

27 

2.4 

2.496 

11 

28 

3.6 

3.713 

11 

29 

48 

4.629 

30 

6.0 

6.213 

31 

7.2 

7.046 

32 

8.4 

8.579 

Information  for  plates 

DAY  TEMP       HA 

15  68  0.108 

1 1  70  0.042 

5  70  0.058 
20  70  23.975 

6  70  0.342 
13  75  0.408 

8  75  0.158 

8  75  0.092 

18  77  0.092 

18  77  0.108 

7  75  0.063 

12  65  0.258 

1 1  65  23.958 
3  60  0.083 

10  65  0.008 

12  60  0.233 

12  60  0.250 
5  60  0.150 
5  60  0.150 
5  60  23.883 

19  60  0.237 

13  68  0.175 
15  68  0.208 
15  68  0.275 

20  70  0.096 
13  68  0.212 
20  70  0.096 
20  70  0.113 
13  75  23.829 

9  75  0.213 

13  75  23.846 

18  77  0.179 


1-32 

Z  D     REFR 

41.15  48.25 
41.14  48.05 
41.14  48.06 
41.14  48.05 
41.21  48.18 
41.26    47.78 

41.16  47.63 
41.14    47.61 

41.14  47.43 

41.15  47.43 
41.14  47.60 
41.18  48.58 
41.14  48.51 
41.14    48.99 

41.14  48.51 

41.17  49.04 

41.18  49.05 

41.15  49.00 
41.15    49.00 

41.15  48.99 

41.18  49.04 

41.16  48.27 

41.17  48.28 

41.19  48.32 
46.14  57.22 
46.16  57.47 
46.14    57.22 

46.14  57.23 

46.15  56.70 

46.16  56.71 
46.15  56.70 
46.15    56.49 
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Table  3:  Refraction  Information  for  plates  33-64 

HA  Z  D  REFR 

0.117  46.14  56.48 

0.167  46.15  57.79 

0.258  46.17  57.82 

0.204  46.16  57.80 

0.329  46.18  58.41 

0.192  46.15  58.35 

0.171  46.15  57.79 

0.088  46.14  58.33 

0.233  46.16  58.36 

0.167  46.15  57.78 

0.146  46.15  57.45 

0.075  46.14157.44 

0.088  51.14  68.46 

0.083  51.14  68.20 

0.250  51.15  68.23 

0.050  51.14  68.20 

23.90  51.14  67.56 

0.354  51.16  67.61 

0.325  51.16  67.61 

0.075  51.14  68.90 

0.308  51.16  68.90 

0.283  51.15  68.89 

0.183  51.15  68.87 

0.250  51.15  69.55 

0.000  51.14  69.52 

0.271  51.15  69.55 

0.154  51.14  69.53 

0.050  51.14  68.46 

0.133  56.14  82.15 

0.350  56.14  80.76 

0.083  56.14  82.62 

0.262  56.14  83.42 


PL 

RAO 

LMST 

MON 

DAY 

TEN 

33 

9.6 

9.717 

1 

18 

77 

34 

10.8 

10.967 

5 

12 

65 

35 

12.0 

12.258 

4 

9 

65 

36 

13.2 

13.404 

4 

12 

65 

37 

14.4 

14.729 

7 

3 

60 

38 

15.6 

15.792 

6 

5 

60 

39 

16.8 

16.971 

4 

19 

65 

40 

18.0 

18.088 

4 

19 

60 

41 

19.2 

19.433 

6 

5 

60 

42 

20.4 

20.567 

10 

11 

65 

43 

21.6 

21.746 

10 

15 

68 

44 

22.8 

22.875 

10 

15 

68 

45 

0.0 

0.088 

10 

18 

68 

46 

1.5 

1.583 

11 

11 

70 

47 

3.0 

3.250 

11 

20 

70 

48 

4.5 

4.550 

1 

7 

70 

49 

6.0 

5.900 

1 

13 

75 

50 

7.5 

7.854 

2 

7 

75 

51 

9.0 

9.325 

2 

7 

75 

52 

10.5 

10.575 

2 

11 

75 

53 

12.0 

12.308 

4 

11 

65 

54 

13.5 

13.783 

4 

12 

65 

55 

15.0 

15.183 

4 

12 

65 

56 

16.5 

16.750 

6 

12 

60 

57 

18.0 

18.000 

6 

12 

60 

58 

19.5 

19.771 

9 

15 

60 

59 

21.0 

21.154 

9 

15 

60 

60 

22.5 

22.550 

10 

18 

68 

61 

0.0 

0.133 

10 

13 

68 

62 

6.0 

6.350 

1 

18 

77 

63 

12.0 

12.083 

4 

12 

65 

64 

18.0 

18.263 

6 

5 

60 

CHAPTER  3 
LEAST  SQUARES 

Since  we  will  use  the  method  of  least  squares  to  determine  our  plate  parame- 
ters for  both  the  single  plate  overlapping  plate  reduction  techniques,  we  will  first 
give  an  outline  of  this  data  analysis  method. 

Theory  of  Least  Squares 

Consider  a  set  of  n  observations  represented  in  vector  form  by  x©,  and  assume 
these  observations  are  unbiased  (i.e.  without  systematic  error).  The  vector  of 
true  values  is  represented  by  x  and  v  =  x  -  xq  is  the  vector  of  errors  in  the 
observations.  From  the  assumptions  that  the  observations  are  unbiased  we  can 
assume  that  the  errors  will  have  a  Gaussian  distribution  (multivariate  normal 
distribution)  and  we  will  regard  their  covariance  matrix,  a,  as  being  known.  Thus 
we  can  express  the  probability  density  distribution  function  of  the  errors  as  the 
multivariate  gaussian 


<^(v)  =  Ce-i^^'^"'\  (37) 


We  further  assume  that  the  observations  satisfy  a  set  of  equations  of  condition 

F(xo  +  v,a)  =  0  (38) 
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where  a  is  a  vector  of  q  parameters.  Some  of  the  components  of  the  vector 
equations  F=0  may  not  explicitly  contain  any  of  the  components  of  x;  such 
equations  would  be  condition  equations  involving  parameters  only.  In  this  analysis 
we  will  distinguish  these  from  the  general  condition  equations  by  calling  them 
parameter  constraints. 

If  the  q  parameters  are  not  mutually  independent,  the  relations  existing  be- 
tween them  must  be  included  among  the  parameter  constraint  equations.  The 
observations  must  be  functionally  independent,  that  is,  not  one  of  the  n  observa- 
tions can  be  derived  from  any  or  all  of  the  remaining  («-l)  observations.  There 
must  be  at  least  n  observations;  in  other  words,  n>q.  Otherwise  there  will  be  a 
deficiency  in  the  model.  If  q  =  n  then  such  a  circumstance  will  lead  to  a  unique 
solution.  When  n  >  q,  redundancy  is  said  to  exist  and  adjustment  is  needed  in 
order  to  obtain  a  the  best  set  of  estimates  for  the  model  variables.  Adjustment  is 
meaningful  only  in  those  cases  in  which  the  data  available  exceed  the  minimum 
necessary  for  a  unique  determination.  Since  the  data  points  are  usually  obtained 
from  observations,  which  contain  errors,  redundant  data  are  usually  inconsistent 
in  the  sense  that  each  sufficient  subset  yields  different  results  from  another  subset. 
No  unique  result  (no  one  vector  a)  will  satisfy  all  the  equations,  but  will  give  rise 
to  (at  least  some  non-zero)  residuals,  v. 

Making  the  basic  assumption,  called  the  principle  of  maximum  likelihood,  we 
assume  that  the  set  of  measurements  which  we  obtain  is  actually  the  most  probable 
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set  of  measurements.  Thus  the  best  estimates  of  the  coefficients  a  are  the  ones  that 
maximize  the  probabiHty  of  obtaining  the  particular  set  of  measurements  which 
we  actually  obtained.  Clearly,  the  way  to  maximize  (p  is  to  minimize  the  value  of 
the  exponent.  The  principle  of  maximum  likelihood  thus  leads  to  the  conclusion 
that  we  should  minimize  v^a~^v. 

The  principle  of  least  squares  consists  in  finding  estimates  of  v  and  a 
which  minimize  the  quadratic  form  v^a~^v  while  simultaneously  satisfying 
F(xo  +  V.  a)  =  0.  Thus  if  the  random  variables  to  which  the  observations  re- 
fer are  normally  distributed,  the  least  squares  method  will  give  identical  results 
to  those  from  the  maximum  likelihood  method. 

Traditional  (Linear  Least  Squares) 

Traditional  least  squares  is  when  we  assume  that  the  observation  errors  are 
of  the  same  precision  and  uncorrected  and  that  exactly  one  observation  occurs  in 
each  equation  of  condition.  Since  the  observations  are  uncorrelated  the  covariance 
matrix  a  becomes  diagonal.  Furthermore,  if  the  observations  are  of  the  same 
precision  the  covariance  matrix  will  be  of  the  form  a  =  (ToI,  where  I  is  the 
identity  matrix.  The  minimum  principle  of  least  squares  becomes  : 

n 

-1„2 


V    (7     -^V  =  2_]  ^     ^Vj' 


':'    n  (39) 

T    -1  ^    V^     2  •    • 

V  a     V  =  —  >    V:  — >  minimum. 
1=1 
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n 

Thus  in  this  case  the  minimum  of  ^  V;-  —  v^v   and     v^<7~'v  coincide.  This 

i=l 

last  case  is  the  oldest  and  the  most  classical,  and  possibly  the  one  that  gave  rise 
to  the  name  "least  squares"  since  in  this  case  we  seek  the  "least"  of  the  sum  of 
the  squares  of  the  residuals. 

In  traditional  least  squares  the  condition  equations  F(x,a)  =  0  usually  have 
the  form 

Aa  +  Xo  =  -V.  (40) 

where  A  =  |^.  In  this  case  the  "adjustment  unknowns"  or  "  adjustment 
parameters",  a,  occur  linearly  in  the  condition  equation,  and  each  equation 
contains  exactly  one  observation.  The  "normal  equations"  for  finding  the  estimates 
of  a  are  easily  found  for  this  case.Multiplying  the  above  equation  by  its  transpose 
gives: 


(a^A^  +  xJ)(Aa  +  Xo)  =  (-vT)(-v) 


a'^A'^  Aa  +  2a'^A^Xo  +  xjxo  =  v'^v 


(41) 

for  matrix  A  and  vectors  a  and  x.  Also,  note  that  v^  v  is  not  a  vector  or  a  matrix, 
but  a  number  and  we  may  therefore  transpose  any  of  the  terms  without  knowing 
its  value.   To  find  the    minimum  of  v^  v,  we  differentiate  with  respect  to  each 


component  of  a  and  set  the  derivative  equal  to  zero 

d(v'^v)       aia'^A'^Aa  +  2aTATxo  +  xJxq) 


=  0 


da  da.  (42) 


2('a'^A'^A  +  xjA)  =0 
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yielding 

a'^A'^A  =  -xJa.  (43) 

After  transposing  and  rearranging,  we  get  the  traditional  least  squares  solution 
in  its  simplest  form: 

a  =  -(A^A)-Ia'^Xo  (44) 

the  solution  of  the  normal  equations. 

Nonlinear  Least  Squares 

The  conditional  and  the  constraint  equations  involved  in  an  adjustment  prob- 
lem can,  in  general,  be  nonlinear.  However,  least  squares  adjustments  are  gener- 
ally performed  with  linear  functions,  since  it  is  rather  difficult  and  often  imprac- 
tical, at  least  at  present,  to  seek  a  least  squares  solution  of  nonlinear  equations. 
Consequently,  whenever  the  equations  in  the  model  are  originally  nonlinear,  some 
means  of  linearization  must  be  performed  on  the  equations.  Series  expansions,  a 
Taylor's  series  in  particular,  are  often  used  for  this  purpose,  where  only  the  zero 
and  first-order  terms  are  used  and  all  other  higher-order  terms  are  neglected. 

When  applying  a  series  expansion,  an  initial  set  of  approximate  values  for  the 
unknowns  in  the  equations  must  be  chosen.  The  choice  of  those  initial  approx- 
imations is  an  important  aspect  of  solving  the  problem  at  hand.  Unfortunately 
there  is  no  concrete  and  unique  way  of  choosing  approximations  that  can  be  ap- 
plied to  all  adjustment  problems.  In  all  cases  an  attempt  should  be  made  to  obtain 
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the  closest  approximations  that  can  be  obtained  by  using  relatively  simple  and 
uninvolved  techniques. 

We  will  now  use  the  more  generalized  least  squares  treatment  to  drop  some 
of  the  restrictions  found  in  the  traditional  least  squares  method.  In  the  following 
development 

1.  Observations  may  be  correlated. 

2.  More  than  one  observation  may  appear  in  an  equation  of  condition. 

3.  The  equations  of  condition  are  assumed  nonlinear  in  both  parameters  and 
observations. 

4.  Nonlinear  constraints  among  the  parameters  are  allowed. 

As  before  the  vector  Xo  is  the  vector  of  observations,  x  is  the  vector  of  true 
values  x=Xo  +  v  ,  where  v  is  the  vector  of  corrections.  In  order  to  find  estimates 
of  the  vector  of  parameters  a  and  x  in  the  generalized  least  squares  we  will  use  the 
method  of  Laplace  multipliers.  At  the  solution  the  numerical  value  of  v^a~^v  = 

5.  These  conditions  are  the  same  as  5*  =  v^a~^v  —  2F^(a,  Xq  +  v)A,  where 
A  is  a  Lagrange  multiplier.  This  means  that  the  values  of  S  and  S*  reach  their 
minimum  at  the  same  values  of  a  and  v. 

As  before,  a  has  q  components,  v  has  n  components  and  F  has  p  components. 
Thus  there  are  only  p+n-q  free  components.   Differentiating  S*  with  respect  to 
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the  components  of  a  and  v  gives 

ds* 


where 


■■pxn 


OF 


d\- 

ds*  T 

-—  =  -2A^A 


dvi        d\2 

dh      Ml 

dv\        dv2 


(45) 


x=xo,a=ao 


Apxq  -    [    g^ 


x=xo,a=ao 


9vi         dv2 


3ai  diL.2 
di2  dh 
3ai       3a2 


afp     dtp 


3Vn 


9L 


(46) 


9v„  -I  i=xo,a=ao 


1^ 
aa.n 

d{2 

5a^ 


5a<, 


(47) 


x=xo,a=ao 


The  following  equations  must,  therefore,  be  rigorously  satisfied  at  the  solution 


.-1, 


a     V  =  X   A       n  equations 
A   A  =  0  q  equations 


(48) 


F  =  0  p  equations. 

This  means  that  A  and  X  must  be  evaluated  at  the  solution  for  a  and  for  x.  When 

one  considers  that  A  has  p  unicnowns,  then  one  can  see  that  there  are  exactly  the 

same  number  of  unknowns  as  there  are  equations. 

The  equations  are  in  general  nonlinear  and  therefore  must  be  solved  by 

successive  approximations.  A  general  rigorous  and  non-iterative  algorithm  for  the 


52 
solution  exists  only  for  the  case  that  the  elements  of  x  and  a  occur  linearly  in  the 
function  F.  When  F  is  nonlinear  in  the  elements  of  either  x  or  a,  or  both,  equations 
which  are  practically  equivalent  to  the  F  and  which  are  linear  in  the  pertinent 
variables  can  be  derived,  e.g.,  in  the  following  way:  assume  that  an  approximate 
solution  of  the  above  equations  has  been  obtained  somehow,  linearize  these 
equations  about  the  approximate  solution  and  by  iterations  determine  corrections 
to  it.  Start  by  taking  an  approximation  ao  for  a  and  using  the  observations  Xq 
as  approximations  for  x.  If  we  write  a  =  ao  +  a,  with  a  being  small  corrections 
to  ao,  the  equation  F=0  can  be  developed  as  a  Taylor  series 

F(x,a)  =  F(x„,a„)  +  (X  -  x,)^?:^^  +  (a  -  a„)^!l^  +  0(2)=  0 

ax  da 

F(x,a)  =  F(xo,ao)  +  Xov  +  Aoa  +  0(2)=0 

(49) 

where  X©  and  Aq  are  evaluated  at  xq  and  ao;  and  also,  let  F(xo,  ao)  =  Fq. 

From  equation  (48)  we  have  v  =  crX^A;  so,  by  substituting  into  the  above 
equations  we  get 


rT 


Fo  +  XoaXo' A  +  Aoq  =  0. 


(50) 


Combining  this  with  A^A  =  0  yields  the  Generalized  Normal  equations 

XoaXj     Ao 
AT  0 


A 

a 

+ 

Fo' 
0 

= 

0 
0 

(51) 


Note  that  A  is  needed  for  computing  the  correction  v.  If  XcrX^  is  nonsingular, 
we  can  eliminate  A  from  the  equation  by  multiplying  the  equation  (50)  by 
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A^(XocrXo)       yielding 

A  J  (XoaX^)     Fo  +  aJa  +  A^(Xo<7Xj  r'^Aoa  =  0.  (52) 

Since,  A^A  =  0,  rearranging  gives  us 

A'^iXoaX^r'Aa  =  -A^{XoaX^r'Fo 

(53) 
a  =  -[A^{XoaXir'A]-'A'^{XoaXlr'Fo. 

Pre-multiplying  equation  (50)  by  (Xoo-Xq  )      and  solving  for  A  gives 

(54) 


{XoaX^r'Fo  +  A  +  (Xo<TXj)-iAa  =  0 


A  =  -(Xo^Xj)-iAa  -  {XoaXlr'Fo. 

A  =  -W(Aq  +  Fo)  (55) 

where  W  =  (Xo<tXo  )     .  W  is  frequently  called  the  weight  matrix. 

Since  we  know  A^A  =  0,  we  can  pre-multiply  equation  (55)  by  A^,  and  get 


A'^A  =  -  A'^WAa  -  A^WFo  =  0.  (56) 


Solving  for  a  gives 

1-1 


Q  =  — 


a'^wa 


A^WFo  (57) 


It  can  be  shown  (Brown  1955;  Jefferys  1980,  1981)  that  (A'^WA)"^  is  the 
covariance  matrix  of  a,  so  that  the  square  roots  of  its  diagonal  terms  are  the 
standard  deviations  of  the  corresponding  terms  of  a.  Note  that  XaX^  is  not 
always  nonsingular;  for  instance,  when  F=0  contains  equations  in  which  none  of 


54 

the  components  of  x  occur  explicitly  i.e.,  parameter  constraints.  This,  however, 
produces  in  general  no  essential  singularities. 

It  is  now  possible  to  solve  for  v.   From  v  =  aX^A  we  get 

v  =  <tX^[WA[A^WA]-^A'^W-W}Fo.  (58) 

Then 

(59) 

"new  =  Xq  +  V 

define  the  improved  solution.  If  the  improved  solution  is  insufficiently  accurate, 
the  process  should  be  repeated  until  satisfactory  convergence  is  attained. 

The  natural  starting  point  for  this  scheme  is  to  set  Xo=  x  as  the  initial 
approximation,  so  that  v  =  0,  and  to  use  a  vector  a  of  initial  approximations 
for  a  for  ao-  Obtaining  a  suitable  a  is  sometimes  a  difficult  question,  it  must  be 
close  enough  to  the  solution  that  the  process  given  above  converges. 

Comparison  with  Traditional  Least  Squares 

It  is  informative  to  compare  this  general  method  with  the  results  from  the 
traditional  least  squares  method.  In  the  traditional  case  we  are  dealing  with  the  a) 
noncorrelated  observations  of  equal  precision  (i.e.  a  =  aol^a  =  (Tolm)  of  which 
b)  exactly  one  observation  occurs  in  each  equation  of  condition,  and  that  c)  the 
equation  of  condition  is  solved  with  respect  to  it,  we  have  X  =  — Im  and  therefore 
[XaX^]      =  ^Im.  so  that  in  this  very  specialized  case  a  =  -(A^A)~  A^Fq 
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Also,  the  observation  errors  reduce  using  this  simpHfication  and 

X  =  -I., 


a  =  Gol 


m 


W=— I, 


(60) 


to 


V  =  <7X^[WA[A^  WA]-iA^  W  -  W}Fo 


V  =  ^oIm(-Im)[— ImA(A^— ImA)-lA^— Im  -  ^U]Fo  (^i) 


<70 


1 
<70 


.(a^a) 


-1 


T 


0- 


This  formula  is  well  known  from  classical  least-squares  theory  and  evidently 
independent  of  Qq. 


CHAPTER  4 
METHOD  OF  PLATE  REDUCTIONS 

We  will  now  use  the  theory  of  least  squares  for  the  determination  of  the 
plate  parameters.  Recall  that  from  chapter  2  the  determination  of  the  equatorial 
coordinates  is  a  two-step  process.  First  the  reference  stars'  equatorial  coordinates 
are  converted  to  standard  coordinates.  Second,  after  the  functional  form  of  a 
plate  model  is  established  between  the  measured  coordinates  and  the  standard 
coordinates  a  least  squares  reduction  is  performed  on  the  measured  and  standard 
coordinates  of  the  reference  stars  to  determine  the  plate  parameters.  Usually  a 
great  deal  of  analysis  is  involved  in  determining  the  functional  form  of  the  plate 
model.  First  a  simple  4  or  6-constant  model  is  used  and  approximate  values 
are  obtained  for  some  of  the  plate  parameters.  A  least  squares  reduction  using 
only  the  reference  stars  on  the  plate  (single  plate  reduction)  is  done  to  find  the 
best  values  for  the  parameters.  After  the  reduction,  the  residuals  from  the  actual 
measurements  (x,  y)  and  those  calculated  using  these  plate  parameters  and  the 
cataloged  equatorial  coordinates  are  studied.  From  the  residual  plots  it  is  decided 
which  new  terms,  if  any,  need  to  be  added  to  the  plate  model.  Once  the  "best" 
plate  model  has  been  established,  the  plate  parameters  corresponding  to  this  model 
are  used  to  determine  the  equatorial  coordinates  of  all  the  stars  imaged  on  the 
plate  (not  just  the  reference  stars).  After  position  estimates  for  all  the  stars  have 
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been  determined  an  overlap  solution  can  be  performed.  This  is  because  in  an 
overlap  solution  all  the  images  of  stars  that  occur  on  at  least  two  plates  or  are 
reference  stars  are  considered  in  the  normal  equations. 

Single  Plate  Reduction:  Using  Least  Squares  to  Determine  the  Plate  Parameters. 

We  will  first  reduce  each  plate  separately.  In  the  next  section  we  will  develop 
a  multiple  plate  reduction  technique.  We  will  use  the  formulation  developed  by 
Eichhom  ( 1985)  to  set  up  the  equation  for  a  single  reduction  and  also  an  overlap 
reduction. 

In  a  least  squares  adjustment  we  regard  the  measured  relative  coordinates 
of  the  images  and  the  position  estimates  from  the  reference  star  catalogue  as 
observables.  These  two  types  of  observables  give  rise  to  two  types  of  condition 
equations.  Consider  data  from  only  one  plate  on  which  there  are  n  stars,  m  of 
which  are  reference  stars.  Without  restricting  generality,  we  may  assume  that  the 
measured  image  coordinates  of  the  i/^^^  star,  x^,  y^  are  related  to  their  standard 
coordinates  ^^i,  //^  by  the  equation 

'''')=s(^^)+Z,a.  (62) 

Si.  is  the  model  matrix,  which  depends  on  ^,  rj,  and  possibly  other  data  pertaining 
to  the  star  in  question,  such  as  magnitude  or  color  index,  a  is  the  vector  of  the 
plate  parameters,  which  constitutes  a  subset  of  the  adjustment  unknowns  and  is 
a  vector  of  dimension  p.    In  what  follows,  we  shall  keep  the  model  matrix  S 
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general  and  hold  to  the  assumption  that  Sa  is  at  least  one  order  of  magnitude 
smaller  than  s^  and  st). 

Using  this  relation  we  establish  the  first  set  of  condition  equations 


IC: 


{■2mxi) 


\\ym 


,772(02,^2) 


--la 


-  i=,,a 


i^ma 


=  0. 


(63) 


In  our  approach,  the  observables  are  x^^  and  y^  and  we  regard  the  adjustment 
parameters  to  be  not  only  the  plate  constants  a,  but  also  the  spherical  coordinates 
(o!fi,6fi)  of  the  stars.  ^  and  r]  are  used  only  as  auxiliary  quantities.  We  consider 
the  tangential  point,  oq  and  6q,  as  being  known.  As  stated  in  Chapter  2  if  the 
assumed  tangential  point  is  not  the  actual  one,  then  this  is  can  be  accounted  for 
by  an  appropriate  form  of  2  and  additional  components  of  a. 

The  second  type  of  observations  (estimates  of  the  reference  star's  spherical 
coordinates  from  an  existing  catalogue)  states  that  the  actual  spherical  coordinates 
(ci^,  Sfj^)  are  equal  to  those  estimated  in  the  catalogue  plus  corrections.  To  take 
care  of  convergence  of  the  meridians,  we  introduce  the  factor  cosS  in  the  right 


ascension  equations  yielding 
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/(Qic  -  qJcOS^i         \ 

(Qjc  ~  ^2)  COS  ^2 
^2C  ~  ^2 


(64) 


(Oimc  -  Qm)  COS  Sm 

where  aic-.S^c  are  the  catalogue  positions  and  ai...dm  are  the  adjustment  un- 
knowns which  occur  both  in  F=0  and  G  =0.  We  will  now  combine  the  4m 
condition  equations  F  and  G  into  one  set  H. 

H'^=  (F'r,G'r).  (65) 

Each  equation,  regardless  of  whether  it  belongs  to  F=0  or  G=0  ,  contains  exactly 
one  observable,  therefore 


S'— 


4m- 


(66) 


We  assume  that  the  measurements  of  all  the  coordinates  are  free  from  bias 
and  not  correlated.  The  variance  (in  units  of  squares  of  the  units  in  which  x^ 
and  y^  were  measured)  of  \^i  is  Ufj,  and  y^  is  ^^.  The  covariance  matrix  of  the 
rectangular  coordinate  measurements  is 


(Tx  =  diag(z/i,  01,  1^2,  (/)2,  ...,  //m,  <^m) 


(67) 


and  in  most  cases,  it  will  suffice  to  write 


O'mL'2 


'X  —  ^m^2m 


(68) 
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where  cr,„  is  the  common  variance  of  all  coordinate  measurements  and  Ijm  is,  of 
course,  the  identity  matrix  of  order  2m. 

Also  let's  assume  that  all  the  reference  star  coordinate's  are  uncorrelated, 
so  that  the  covariance  matrix  aai>  of  the  catalogues  reference  star  coordinate 
estimates  qcos(5  and  6  will  be 

^at/ =  diag(^i,d-i,...,^„,<Tm).  (69) 

The  factors  cos^Si/   are  introduced  in  order  to  make  s/p^  the  dispersion  of 
auCOsiSiy),  that  is  ,  the  dispersion  on  the  great  circle.  Not  only  will  this  simplify 
some  of  the  formulas  we  need  in  our  calculation,  but  it  compensates  for  the 
convergence  of  the  meridians. 
Thus,  the  observations  are 


T' 


and  their  covariance  matrix  is 


cr  =  diag(i/i,0i,z/2,(?i)2,...i/^,(^^;Pi,cri,p2.cr2,.../9^,(T,„)  (71) 


or 

Vx      0 

0      a. 


P^      0 

0      a^ 
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The  vector  a  of  adjustment  parameters  consists  of  the  corrections  to  the 
reference  star  positions  and  the  initial  plate  parameters.  As  approximations  for 
the  position  we  choose  their  cataloged  values  which  renders  Go  =  0.  As  initial 

approximations  to  the  vector  of  the  plate  parameters  we  choose  0.  Thus  we  have 

Fo  =FiacJc,0]=d 
/Fo 
Ho  = 

\Go/        \0 


(73) 


The  vector  of  the  adjustment  parameters  a  =  | ,  with 

/?    =  (cos(5idai,d^i cos6mdam,dSm)-  (74) 


From  our  least  squares  analysis  in  Chapter  3,  the  corrections  to  the  parameters 
Q  are  given  by 

Q  =  -  [  A"^  WA]  -  ^  A^  WHo  (75) 

where  W  =  [XaXJ)         Since  all  the  equations  have  only  one  observation, 
X  =  I  and  our  weight  matrix  becomes  W  =  cr"^ Im  so  we  can  write 

Q  =  ~[A^a-^A]-^A^a-'^Ho.  (76) 

In  order  to  solve  the  normal  equations,  we  have  to  evaluate  explicitly  the 

(■'\ 

terms  of  the  matrix  A  and  since  a  =  I       I ,  we  can  rewrite  the  above  equation  as 
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A  is  the  Jacobian  matrix  of  the  condition  equations  with  respect  to  the  adjustment 

parameters,  so  we  can  generally  write, 

' dH\  _  (  rXF^\  _  (  {^)      (f) 


If  |q|  <<  1  and  if  the  initial  parameter  vector  is  a  =  0  then,  the  values  we 
obtain  for  a  from  equation  (77)  will  be  the  actual  parameter  values  rather  than 
the  corrections. 

Let's  take  a  closer  look  at  what  these  partial  derivatives  equal; 


.9l^nJ  \d{aft  cos 6^,  6^)  J  (79) 

^here  B  =  (||^)  (^^  ^^^  When  u  does  not  equal  ^,  (^)  =  0. 
Therefore  (|j)  is  block  diagonal,  consisting  of  m  blocks  of  dimension  2  x  2  on 
the  main  diagonal;  these  blocks  are  -sBi/.  The  factor  matrix  is  required  because 
we  regard  cos6  da,  and  not  da  as  the  correction,  thus 

The  other  elements  of  A  are 


dF 
da 


(81) 


V-sJ 


(dG,  dG,  dG,  dG,  . 

cosT75a7  dSi  ■  ■  ■       cos  bnOan  d6n 

dG  dG^  dG.  dG^    I  (^2) 

C0s6i^ai  "SFT  '  '  '         cos  dnOOln  ddn    ' 
—   Aon 
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and 


da 
Therefore, 


dG 

0.  (83) 


A  =  -(f     ;)  (84) 

We  will  substitute  this  into  equation  (77)  and  examine  the  remaining  terms. 


A^a-^  = 


.sB^     I  \  /aji       0 


X 


E^a-i        0 


-1 


(85) 


r'BTa-iB  +  a-1     .BV'SA  ^^^^ 

It  is  informative  to  look  at  the  dimensions  of  A  and  the  elements  that  make  up 
A.  For  m  reference  stars  and  p  plate  parameters  we  have 

•  B:  block  diagonal  (2m  x  2m),  whose  elementary  matrix  is  a  (2  x  2). 

•  <^^^  '■  block  diagonal  (2m  x  2m),  whose  elementary  matrix  is  a  (2  x  2  ). 

•  cr^^  :  block  diagonal  i2m  x  2m),  whose  elementary  matrix  is  a  (2  x  2). 

•  S  :    {2m  X  p) 

Therefore  the  dimensions  of  the  2  x  2  blocks  of  A  are: 

•  (1.1)  position  :   {2m  x  2m) 

•  (1,2)  position  :    {2m  x  p) 

•  (2,1)  position  :    (m  x  2p) 


64 

•      (2,2)  position  :    (am  x  m) 

After  substituting  in  equation  (86)  into  equation  (77)  we  get 

and  expanding  this  gives  two  equation  vectors 

(.s-B^cr-^B  +  a-^)^  +  sB'^a-^Ea  =  .sB^^a^M  (88) 

and 

sE'^(7^^B3  +  2^(7~'Sa  =  E^a'^d.  (89) 

Solving  the  first  for  3  gives 

/?  =  (s-B^fT^ifi  +  ir^^y  sB^^^d  -  Sa).  (90) 

Substituting  in  ,i  to  the  second  equation  and  letting 

Y  =s-a-^B(^s-B^a-^B  +  a'^Y^B'^a^^  (91) 

we  (after  some  algebra)  get  the  following 

hT{Y  -  a-')d  =  sT(Y-a-i)Sa.  (92) 


Let  J  =  cr~^  -Y  or 


J  =  (t-1-5-V^^b(3-bT(j-1B  +  a-^]     B'^a 


-1 

X 


(93) 
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using  the  inversion  lemma  (Henderson  and  Searle,  1981)  we  get 

J=  [a^  +  s-Ba^B'^^      .  (94) 

Substituting  this  equation  for  J  into  equation  (92)  yielding 

E^J2a  =  2^Jd.  (95) 

We  can  solve  this  equation  for  a: 

a=  (z^JSy  E^Jd.  (96) 

Note  that  J  is  also  block  diagonal,  consisting  of  2  x  2  blocks  on  the  main 
diagonal.   We  can  therefore  write  more  explicitly 

a  =  [  ^  Sjj.E^  j      Yl  -iJ^J'^d^.  (97) 

Next,  we  have  to  find  explicit  expression  for  the  terms  of  J.  These  depend  on  o^^ 
Qa  and  the  terms  of  B.  The  terms  of  B  are  essentially  the  terms  of  (  qT^^X  )  ■  Since 
we  have  good  approximations  of  the  values  of  i,  rj  which  are  approximately  {\/s) 
and  (y/s)  respectively,  we  will  determine  J  from  these  quantities.  Following  the 

procedure  suggested  by  Eichhorn  (1985)  we  will  define 

r2  =  1  +  (''  +  rj^ 

T    =  (,' +  {cos  So  —  r]  s\n  6o)' 

S  =  ^R-sin<^o/T 

(98) 
U  =  (^R(sin  ^0  +  ^  cos  ^o)/T 

V  =  R''(cos<5o  —  r]  sinio)/T 

VV  =  ^^Rcos(5o/T  +  V/R 
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and  note  that 


R'  =  SU  +  VW.  (99) 


B  can  now  be  written  as 


B=(u     ;')  (100) 

From  our  earlier  discussion  we  wrote  J  =  {a^  +  s-Ba^B^)  .  However  witii  the 
formulas  above  ,  the  elements  of  B  can  be  easily  calculated  without  the  explicit 
knowledge  of  the  star's  spherical  coordinates.  If  we  further  introduce  the  variables, 

QS  =  -joUW  +  crSV 

Y  =  p\]-  +  aV^  (101) 

Z  =  pW-  +  aS^ 
and  introduce 


s:)=---=(v  Xi 


the  weight  matrix  J  can  be  written  in  the  form 

Note  that  for  any  nonsingular  matrix  D  the  following  is  true 
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Therefore  the  determinant  of  J  is  A  =  i/o  +  Zo  +  Yu  +  R^ap  and  we  get  the 
weight  matrix  belonging  to  the  z/'*^  star  as 

j^  = \ f^'^  +  Y,        Q,     \ 

We  now  need  to  explore  more  fully  the  computation  for  the  corrections  to 
the  star  positions.  We  now  assume  that  we  have  calculated  the  vector  of  plate 
parameters,  a,  from  equation  (97),  and  henceforth  we  consider  that  the  vector 
d  —  Ea  is  known.  As  mentioned  before,  we  have  m  different  independent  sets 
of  equations  for  calculating  the  corrections  /^J  —  {cos 6 ftda^^dS^)  for  each  star. 
From  equation  (90)  we  have  the  formula  for  the  corrections  to  the  star  parameters, 
expanding  this  equations  out  (and  substituting  a~^  =  s'^^a^^  we  have 

8  =  -(b'^^-^B  +  a-^y\^B'^a-^id  -  Sa) 

^  I DK-  m  i)c^  n^ 

(106) 
Remember  that  3  is  different  for  each  star  and  that  the  symbols  W,  U,  V,  S 
above  are  simple  numbers  and  not  matrices.  To  calculate  the  inverse  of  the  2  x 
2  matrix  in  this  equation,  we  need  its  determinant.  After  some  algebra  we  get 
the  determinant  of  3  to  be 

A  =  — — {i^<p  +  Yu  +  Z<f>  +  R^piy).  (107) 

i/opcr  ' 
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Except  for  the  factor,  ^j^t^,  this  is  the  reciprocal  of  the  determinant  of  J.  After 
multiplying  the  equation  for  J  (equation  106)  by  its  determinant  we  get 
1 


3  = 


.s{vo  +  Zo+  Vv  +  R^pa) 
P     0 


0     a 


% ;;  :  D-^t :  :•  ^- 


(d-Sa). 

(108) 


Note  that  3  is  in  radians.  Its  is  (if  one  investigates  more  closely  )  a  sophisticated 
weighted  average  of  the  catalogue  positions  and  that  which  follows  from  the  plates. 

Final  Positions  from  a  Single  Plate  Reduction, 

Once  we  obtain,  by  an  adjustment  procedure,  estimates  for  the  plate  constants, 
we  use  them  with  an  inverse  of  the  plate  model  to  determine  estimates  for  the 
standard  coordinates  for  all  the  measured  images  (reference  and  non-reference). 
Then  using  the  inverse  gnomonic  projection,  the  equatorial  coordinates  are  de- 
termined. One  of  the  drawback  of  the  single  plate  solution  is  that  if  the  star 
occurred  on  more  than  one  plate,  the  estimates  for  the  same  star's  coordinates 
will,  of  course,  be  different  depending  on  which  plate  they  originate  from.  Thus, 
the  single  plate  solution  does  not  use  the  simple  fact  that  a  star  can  at  one  time  oc- 
cupy only  one  position.  Since  all  the  plates  are  reduced  individually,  a  star  that  is 
found  on  n  plates  will  have  n  different  values  for  its  coordinates.  Clearly  the  star 
only  has  one  position  at  one  time,  but  this  cannot  be  obtained  without  mathemat- 
ically enforcing  this  constraint  on  a  multiple  plate  reduction  process.  Instead,  the 


69 

separately  computed  values  are  averaged  to  give  the  "best"  (in  the  mathematical 
statistics  sense)  estimate  for  the  stellar  position.  In  addition,  the  random  errors  are 
reduced  by  a  factor  of  -4^  but  this  does  not  reduce  the  systematic  errors.  Another 
drawback  of  the  single  plate  solution  is  that  when  using  only  the  reference  stars 
in  the  reduction  we  would  more  than  likely  be  extrapolating  the  magnitude  terms 
to  the  field  stars  which  are  typically  fainter  than  the  reference  stars. 

It  was  considerations  of  this  kind  which  led  to  the  method  of  overlapping 
plates,  in  its  full  form  first  published  by  Eichhom  (1960).  The  method  of 
overlapping  plates  avoids  the  multiplicity  of  several  "best"  star  positions  and 
increases  considerably  the  accuracy  of  at  least  some  of  the  plate  constants,  thus 
reducing  the  systematic  errors  (i.e.,  increasing  the  accuracy)  of  the  final  positions. 
However,  its  implementation  requires  an  arithmetic  effort  exceeding  that  which 
is  necessary  for  traditional  plate  reductions  by  a  factor  of  at  least  2  orders  of 
magnitude.  This  task  has  become  possible  only  within  the  last  few  years,  because 
cheap,  powerful  computers  have  become  easily  available,  removing  the  label  of 
"major  job"  from  a  complete  overlap  solution. 

Overlap  Plate  Reduction 

The  basic  principle  of  the  overlapping  plate  method  is  that  all  the  star 
coordinates,  not  just  the  reference  stars,  are  considered  together  with  the  plate 
parameters,  as  adjustment  parameters  (unknowns).    Since  the  overlap  solution 


70 

regards  the  coordinates  of  each  star  as  adjustment  parameters  as  long  as  images  of 
the  star  were  measured  on  more  than  one  plate,  the  number  of  parameters  increases 
in  proportion  to  the  ratio  of  the  field  stars  to  the  reference  stars.  Obviously  — 
since  the  number  of  unknowns  has  increased  dramatically  —  the  analytical  and 
computational  problem  is  exponentially  more  complex  and  time  consuming  than 
the  single  plate  reduction.  In  addition,  the  adjustments  pertaining  to  individual 
plates  also  lose  their  mutual  independence  because  any  two  plates  which  have  a 
star  in  common  can  now  no  longer  be  reduced  individually. 

Another  unavoidable  defect  of  the  single  plate  reduction  is  a  large  parameter 
variance  (Eichhom  and  Williams,  1963).  Even  if  the  geometry  of  the  imaging  sys- 
tem was  known  precisely  and  the  errors  of  the  measuring  machine  were  modeled 
exactly,  there  will  still  exist  unavoidable  random  errors  in  the  x-y  measurements 
and  reference  star  positions.  This  leads  to  an  unfavorable  error  propagation  of 
the  resulting  star  positions,  in  particular  when  the  number  of  parameters  is  large 
and  the  number  of  reference  stars  per  plate  is  small.  The  accidental  errors  of 
these  plate  parameter  estimates  thus  reappear  as  systematic  errors  in  the  positions 
that  were  calculated  with  them.  Such  systematic  differences  even  showing  up 
as  functions  of  model  terms  are  unavoidable  (and  their  magnitudes  can  even  be 
estimated)  for  those  stars  (field  stars)  which  were  not  involved  in  the  derivation  of 
the  single-plate  parameter  estimates.  An  overlap  solution  enforces  the  constraint 
that  a  star  can  only  have  one  position  at  one  time,  producing  more  accurate  plate 
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parameters,  thus  inevitably  leading  to  smaller  systematic  errors  in  the  star  posi- 
tions. Also,  since  the  positions  themselves  are  regarded  as  adjustment  unknowns, 
there  can  not  be  any  systematic  residuals  between  plates  (except  those  introduced 
by  model  deficiencies).  In  addition,  when  there  is  a  strong  overlap  between  plates, 
fewer  reference  stars  are  needed  per  square  degree  to  obtain  the  same  accuracy 
since,  in  effect,  the  total  area  covered  by  the  overlapping  plates  is  treated  as  one 
large  plate. 

We  assume  that  there  are  altogether  m  stars  involved  in  the  adjustment, 
including  both  reference  and  field  stars  found  on  more  than  one  plate.    Each 

star  number:  /i  =  1,  2 m.  The  numbers  ji  are  assigned  to  the  stars  in  some 

organized  fashion.  There  are  n  plates  and  each  plate  is  assigned  a  number  i/  = 
1,  2,...,  n.  On  the  v^^  plate  there  are  mu  stars.  For  example,  (x^j/,  y^^y)  are  the 
measured  rectangular  coordinates  of  the  images  of  the  /x'*^  star  on  the  u^^  plate. 
mr  of  the  m  stars  are  reference  stars  and  the  reference  stars  are  assigned  a  number 

f^r  11  I-I-T21  ■■■1  fJ^Tmr  • 

The  stars  used  in  a  overlap  reduction  are  either  field  stars  which  occur  on  at 
least  two  plates  or  reference  stars.  The  estimates  for  the  field  stars'  equatorial 
coordinates  are  found  from  the  single  plate  solution.  When  the  star  appears  on 
more  than  one  plate  the  average  of  the  single  plate  solution  is  used.  Thus  the 
averaged  a^,  6^  of  the  field  stars  are  used  to  find  the  standard  coordinates  <^^i/, 
Tji^ii/  on  all  plates  on  which  the  star  appears.    However,  for  the  reference  stars 
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the  catalogued  values  Qc/^,  Scf_i  are  used  as  initial  approximations  to  these  stars' 
spherical  coordinates. 

For  the  establishment  of  the  condition  equations,  it  is  important  to  fix  the 
order  of  the  equations  and  the  order  of  the  adjustment  parameters.  The  condition 
equations  are  arranged  in  the  following  order.  Plate  by  plate  we  establish  the 
F-type  equations  for  each  star  on  that  plate  which  either  occurs  on  at  least  on 
other  plate  or  as  a  reference  star  in  order  of  increasing  star  number.  These  are 
followed  by  the  G-type  equations  which  are  produced  by  the  estimates  of  the 
spherical  coordinates  of  the  reference  stars,  also  in  the  order  of  increasing  star 
numbers.  The  order  of  the  observations  is  the  same  as  the  order  of  the  equations 
which  they  generate.  We  chose  the  order  of  the  parameters  as  follows:  star  by 
star,  in  numerical  order,  corrections  to  the  coordinates,  i.e.  cos^^c/q^  and  dS^j, 
with  fi  =  \,  2,...  m  are  followed  by  the  frame  parameters  ai,  ...,  an. 

The  equations  of  condition  are  identical  to  those  in  the  single  plate  solution. 
As  before  the  total  vector  of  conditions  equations  H=0  consists  of  "plate  equa- 
tions", F=0  and  "catalogue  equations",  G=0.  First  let's  discuss  the  F=0  equations. 
Plate  by  plate  we  establish  the  F-type  equations  for  each  star  on  that  plate  which 
either  occurs  on  at  least  one  other  plate  or  as  a  reference  star.  The  observations 
are  x^j^and  y^j,,  the  measured  coordinates  of  the  image  of  the  fx^^  star  on  the  j/"^ 
plate.   Exactly  one  observation  occurs  in  each  equation  of  condition,  as  before. 
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The  frame  equations  for  the  u  plate  are 


F.  = 


^/jf 


'  I-IU 


t,ui/ 


.'Irn^ 


—  uu^u- 


(109) 


There  are  my  such  equations  for  each  plate,  for  a  total  of  Irrii,  equations,  iau 
and  "qui,  are  the  standard  coordinates  of  the  /f"^  star  with  respect  to  the  assumed 
tangential  point  of  the  t/"^  plate.  2,,^  is  the  model  matrix  for  the  //'^  star  on  the 
u  "^  plate  and  a^/  is  the  vector  of  plate  parameters  on  the  u^^  plate.  Expanding 
the  above  equation,  we  get 


^{l^^V^V 

—   ^Hu^I^^Li 

-m^"    ~  ^ii'-^'m^^   ^ 

en 

~   '^^t-m^'^^" 

(110) 


At  first  glance  the  above  notation  looks  rather  confusing.  However,  x^^^  ^  is  the 
first  star  on  the  v  '^  plate.  Since  all  the  stars  are  given  a  particular  number 
designation  numbered  1,  2....,  m.  There  has  to  be  a  way  to  keep  up  with  the 
running  star  number  ^  and  the  star  number  on  each  plate,  which  we  will  refer  to 
as  1^1^.  The  number  of  stars  on  the  v^^  plate  run  from  ^^^  to  /u^,^^  These  subscripts 
are  useful  in  designating  the  total  running  star  number  and  the  star  number  on 
a  particular  plate.  See  the  example  at  end  of  this  chapter  on  how  to  assign  the 
running  star  number. 
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The  entire  vector  F  is 


F  = 


Vf„  / 


0. 


(Ill) 


Each  vector  Fi/  has  the  dimension  2miy  x  1  and  the  dimension  of  the  total  vector 

n 

is   Yl  i-^i^  X  !)• 

Now  on  to  the  reference  star  condition  equations.  There  are  a  total  of  rrir 
reference  stars,  numbered  //ri,^r2' •••i/^r„^-  For  each  reference  star  we  have  an 
estimated  value  of  its  position  (Oc^ ,  ^c^t  )  from  the  catalog.  The  G=0  equations, 
look  like  this: 


G 


(acfi,-!  -  ^Mri )  COS  6^^,  +  e,,,, 

(aCMr2    -   a^.2  )  COS  (5^^,    +t^^2 
^CHr2    ~  ''|tr2    +  ^|ir2 

[Occilrmr    "   "/'rm.  )    COS  6 ^i^^^    +  C^, 


^CHrnlr  ^  lirmr       '     ^ 


t^rrrtr 


(112) 


Since  only  a  small  fraction  of  all  stars  are  reference  stars,  it  follows  that 
Hrk  and  tirk+\  of  two  neighboring  reference  stars  normally  do  not  have  the 
same  number.  That  means  that  the  subscript  r  is  just  to  keep  track  of  reference 
stars,  but  just  because  a  star  has  a  subscript  r  it  does  not  mean  that  it  is  given  a  dif- 
ferent number  from  that  of  a  field  star.  G  does  not  consist  of  a  separate  equations 
for  each  plate,  as  they  do  for  the  equations  F,  it  has  a  dimension  of  2mr  x  1 . 
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As  before,  the  total  vector  H  consists  of  the  plate  equations  F  and  the  catalog 
equations  G. 


(113) 

Likewise  the  vector  of  residuals  is 

Hj=(dT,o)  (114) 

The  covariance  matrix  is  similar  to  that  in  a  single  plate  setup,  except  that  cr^  is 
broken  down  into  blocks  belonging  to  the  various  frames.  Again  we  are  assuming 
no  correlation  between  observations,  so  the  covariance  matrix  is  diagonal.  It 
consists  of  v^i/,  (pjiu,  the  variances  of  \nu  and  y^v,  and  the  p^i  ,a^  which  are 
cos^<5^  times  the  variance  of  a^^  and  the  variance  of  S^i,  respectively,  for  those 
/i  which  belong  to  a  reference  star 

C7  =  diag(crxi,(7x2,...,crx„;<7a)  =  f  q"     ^    J.  (115) 

The  dimensions  of  each  (Tx;.  :   'Irui,  x  2m^,  so  the  dimension  of  the  total  matrix 

n 

<^x  '■    ^  (2mi/  X  2m,/),  while  the  dimension  of  cr^  :   'Irur  x  2mr. 
t/=i 

As  with  the  single  plate  each  equation  contains  only  one  observation  so, 
-^  ~  (§x)  =  ^'  wherex  =  (xi,yi,X2,y2,  ■••,Xm,ym;  c>;icCos6i,^i,  ...amcCos(5r„,6m), 
so 


XaX    =  a  = 


0       <7a 

The  mathematical  development  of  an  overlap  solution  is  very  similar  to  the 
single  plate  method  discussed  in  the  last  section.    The  major  difference  is  the 
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increase  in  the  observations  (all  the  stars),  which  enormously  increase,  the  size 
and  structure  of  the  matrices.  We  will  present  the  mathematical  development 
similar  to  the  single  plate  development,  noting  the  changing  the  structure  of  the 
matrices. 

As  before,  the  correction  to  the  parameters  is  given  by 


a  =  -(A^a-'^A]      A^^rj-^Ho 


(116) 


3 


where  q  =  I  1  and  3  is  the  vector  of  corrections  to  the  spherical  coordinates 
(q^  cos 6^1, S^i).  The  total  vector  of  plate  parameters  for  all  the  plates  is  a,  which 
is  made  up  of  the  individual  vectors  for  each  plate 


(117) 


Van/ 
Our  primary  concern  is  the  establishment  of  the  matrix  A  for  the  overlap  case: 


<9H\  _  /^5(F,G) 


(S)    (f) 


daJ-\m.)J-\^^)      (If) 


Similar  to  the  single  plate  reduction  we  have  for  each  star 


(118) 


(119) 
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where  each  B,,^  j,  is  a  2  x  2  matrix.  For  each  plate  we  get 
OF, 


dfi 


—  — o>Bj/ 


0      0      B 


—  s 


r 

0 


Vo 


^l„^^' 


0 


0      B 


Hu^U 


B 


/'^m^  " 


^ 


0/ 


(120) 


where  the  dimension  of  each  Bi/  matrix  is  {2miy  x  2m).    Thus  each  line  pair 

contains  exactly  one  B^^^  matrix  (of  dimension  2x2)  but  only  in  those  column 

pairs  which  belong  to  a  star  whose  image  was  measured  on  the  //'^  plate.  In  other 

words,  the  rows  correspond  to  the  stars  on  plate  //,  while  the  columns  are  for  all 

the  stars,  therefore  B,,^  u  is  not  the  null  matrix  only  when  the  fi^^  star  is  on  the 

v^^  plate.  The  entire  matrix  B  for  all  the  plates  is  written  as 

/sBA 


sB  = 


sB, 


(121) 


V'SBn/ 
n 

and  this  matrix  is  of  dimension  ^  2m^  x  2m. 

1/=! 


Furthermore,  we  have 


5FV 


(122) 


where  2^  is  of  dimension  2miy  x  /jy,  /,y  being  the  number  of  components  of  the 
vector  a  on  the  v^^  plate.  In  Sj,,  each  star  with  number  fi  that  occurs  on  the  u^^ 
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frame  generates  a  2  x  /;/  matrix  S^^.   Thus 


(123) 


^'^H^r,"    ' 


and  the  total  matrix  2  is  of  the  form 

/2. 


da 


Vo 


2  J 


(124) 


and  has  dimension  YL  '~^v  ^  Y^U>- 


Furthermore, 


5G^ 


dO,^       dG, 


=  -K. 


(125) 


d/S        cos  Sfn,da^v  38 ^„ 
All  nonzero  elements  of  the  matrix  K  are  1 .  However,  since  not  all  stars  are 

reference  stars,  K  is  not  the  identity  matrix.  K  has  the  dimension  of  2mrX  2m.  If 

the  //"'star  is  also  a  reference  star  then  the  row  pair  corresponding  to  this  reference 

star  and  the  column  pair  for  this  //'*'  star  will  contain  the  I2  identity  matrix. 

Since  G  is  independent  of  the  plate  parameters,  we  have  (as  in  the  single 

plate  result) 

dG 


da 


=  0. 


(126) 


Combining  all  these  results  (as  in  the  single  plate  reduction),  the  matrix  A  can 
be  written  as 

sB    2 


A-  - 


K      0 


(127) 
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(    n  \ 

and  has  dimension  21  ^  m^  +  nir  1  x  2m  +  ^  /j,. 

Now  tiiat  we  have  looked  at  the  various  matrices  that  make  up  the  matrix  A, 
let's  look  at  the  structure  of  the  normal  equations.  The  structure  of  the  normal 
equations  looks  the  same  as  in  the  single  plate  solution: 

However,  these  equations  are  much  more  complicated  than  in  the  single  plate 
case.    In  detail 


t,-i._/^^bV     KT^rMAB     S 


^   ""  ^=  '  ST,-i  0       A  K      0 


(129) 


which  can  also  be  written  as 

/.-' f  bT<7;1B,+KV'K     .B^ar^H,      ...     sBT<t-1S„\ 
.sH^ariB.  sT'^riBi      ...  0 


A'^<j-1a= 


L  5BT,7r^S„ 


5S'^(T-^B      2'^(T-^S 


(130) 


with  L  =  K^a-^K  +  ^^B^^a-iB. 

Now  let's  take  a  closer  look  at  structure  of  the  matrices  that  make  up  equation 
(130).  The  example  given  later  on  will  help  to  clarify  these  matrices,  so  the 
reader  may  want  to  refer  to  this  example.    Because  ctx  is  a  diagonal  matrix, 
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multiplication  by  this  matrix  leaves  the  structure  of  the  matrix  by  which  it  is 
multiplied  unchanged. 

We  noted  from  the  discussion  above  that  the  Bj/  (2mi/  x  2m)  are  very  sparse 
matrices  and  have  a  structure  similar  to  that  of  K.  Each  row  pair  contains  exactly 
one  non-zero  2x2  matrix,  however  not  every  column  pair  has  one  of  these 
blocks,  (i.e.  some  of  the  column  pairs  contain  all  zero  elements).  Nonzero 
elements  occur  in  a  column  pair  of  Bj/  only  if  the  corresponding  star  occurs  on 
the  plate  corresponding  to  the  row  for  that  plate.  Compact  blocks  with  nonzero 
elements  in  B^jy  will  interact  with  compact  blocks  in  B,y,  wherever  the  /x"^  star's 
image  was  measured  on  the  t^'^  plate.  The  result  of  the  interaction  will  result  in 
a  2  X  2  block  in  the  /i'*^  row  pair  —  and  column  pair  —  of  L.  The  structure  of 
B^(7~^B  is  thus  similar  to  that  of  B^,  except  that  the  u^^  column  pair  would  be 
"replaced"  by  a  matrix  of  the  width  l;y.  The  matrix  K^a~^K  will  be  a  diagonal 
matrix,  however  some  of  the  elements  on  the  diagonal  may  be  zero.  The  result  of 
adding  this  matrix  to  s~B^a~^B  is  that  the  diagonal  blocks  {pn,  a^)  will  add  to 
the  diagonal  blocks  in  the  /^"^  double  row  and  column.  L  is  thus  a  block-diagonal. 

Since  the  dimensions  of  the  matrices  B^j/  are  2m  x  2mi,  and  those  of  Sj, 
are  2miy  x  liy.  the  product  sBjcr^^S^  is  a  sparse  matrix  of  dimensions  2mxlj/, 
whose  nonzero  2x2  elements  sBj,,a~}S^^  occur  only  on  those  of  its  line  pairs 
whose  numbers  correspond  to  a  star  which  occurs  on  the  v^^  plate. 
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S^a'^S  is  block-diagonal,  the  individual  blocks  being  Sj(T~^2t,,  the  matri- 
ces of  the  product  sums  of  each  plate's  model  matrix  with  the  covariance  matrix. 
Finally,  consider  the  right  hand  side  of  the  normal  equations  (equation  128), 

T    -1  " 

A '■a  ^H.  The  resulting  vector  has  the  dimension  (2to  -|-  ^  /j,)  x  1  which  is  a 
vector  of  the  order  of  twice  the  number  of  stars  plus  the  total  number  of  plate 
parameters. 


T_-l 


AV-^H=-  I  Z 


BT     KT\/<t-i       0    \/d 
'T       0    J      0       ar'     [o 


a 


n 

Consider  the  2  x  1  vector:  5  ^  Bj^a'^dfj,^.  Nonzero  contributions  to  these 
sums  are  made  only  when  images  of  the  /.i^^  star  were  measured  on  the  i/^^  plate. 
Substituting  equation  (130)  and  (131)  in  the  normal  equations  (128)  yields 


The  first  row  gives  L/^  +  sB^cr~^Sa=  sB^crjM  and  the  second  row  gives 
6S^(T~^B/i  +  S^cr~'2a  —  H^cr~M.  As  before,  the  solution  to  the  normal  equa- 
tions is  done  in  steps.  First  we  eliminate  the  star  parameters,  /?,  and  solve  for  the 
plate  parameters  a.  Solving  the  first  equation  for  /?  gives  : 

/S^sL-^B'^a-'^id-Sa).  (133) 

Substituting  /3  into  the  second  row  of  the  normal  equations  yields 

2'^(a-i  -.s-cr-^BL-^B'^(7-^)(d-Sa)  =  0.  (134) 
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Let 

J'  =  ^-1  _  .s^a-'BL-iB'^a-i  (135) 

Then 

S'r(J')(d-Sa)  =  0 

STj'd  =  2Tj'2a  (136) 

a=(ETj's)"'E:'rj'd. 
The  inversion  lemma  cannot  be  used  to  simplify  J'  as  it  was  used  to  simplify 
J  because  the  2m  x  2m  matrix  K^cr^^K  term  in  L  is  a  diagonal  matrix  with 
nonzero  elements  on  the  diagonal  only  in  those  double  blocks  whose  ordinal 
numbers  correspond  to  reference  stars,  thus  it  is  singular.  Unlike  the  matrix  J',  J 
is  not  block-diagonal.  Since  we  have  established  above  that  L  is  block-diagonal 
with  blocks  of  dimension  2x2,  the  structure  of  J'  is  the  same  as  that  of  BB^. 
Thus  the  resulting  S^J'2  is  "banded-bordered"  (De  Vegt  and  Ebner,  1972)  and 
thus  there  are  simplified  routines  to  invert  it. 

Now  that  we  have  determined  (i.e.,  estimated)  the  plate  parameters,  we  can 
use  them  to  find  the  star  parameters  by  equation  (133).  In  terms  of  the  individual 
stars 

(n  \  -1    n 

^^,^r,%l  +  ^"  E  bJ.'^;.'B^,         Y1  Bj.^;.'(d^,  -  S^,a,).  (137) 
j/=i  /       i/=i 

The  symbol  6^i,  is  of  Kronecker  type,  it  equals  zero  if  the  ^"'  star  is  not  a  reference 

star,  and  it  equals  1  for  reference  stars.  The  2  x  1  vector  l3fi  =  (do^  cos  6fi,  dS^)^ 
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is  a  sub-matrix  of  the  previously  encountered  2mi/  x  1  vector  Jj/  which  stands 
for  the  pairs  (da  cos  6.  dS)  of  all  the  stars  which  occur  on  the  v^^  frame.  Thus  if 
a  star  occurs  on  n  plates,  there  will  be  n  different  /i's  after  the  reduction,  whereas 
in  the  case  J^,  all  the  star  parameters  for  one  star  are  calculated  after  all  the 
plate  parameters  have  been  calculated.  This  substantially  amounts  to  taking  a 
weighted  mean  for  all  the  contributions  of  all  frames  on  which  images  of  the  star 
were  measured. 


Example  of  Overlap  Reduction. 

Although  the  solution  of  the  normal  equations  in  the  overlap  case  is  very 
similar  tot  that  of  the  single  plate  case,  the  interaction  of  the  stars  imaged  on 
more  than  one  plate  complicates  the  reduction.  It  is  therefore  extremely  useful  to 
illustrate  the  method  with  an  example. 

We  will  use  a  very  simple  example.  Consider  the  case  where  we  have  8  stars 
imaged  on  2  plates,  6  stars  on  the  first  plate  and  6  stars  on  the  second.  Let  three 
of  the  stars  be  reference  stars. 


plate  1 
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2 
* 


plate  2 


5 


4 
* 


6 

• 

* 

• 

:  field  star 

* 

:  reference  star 

Figure  9:  Example  of  two  overlapping  plates 

One  can  see  that  each  plate  has  a  star  in  the  comer  to  which  no  number  has 
been  assigned.  Each  of  these  two  stars  occurs  on  only  one  plate,  and  none  is  a 
reference  star,  thus  in  the  overlap  reduction  they  are  ignored.  After  the  reduction 
has  been  performed,  the  plate  parameters  are  used  to  calculate  the  equatorial 
coordinates  of  these  isolated  stars. 

Using  our  notation  established  earlier  we  have 


m  =  6 
n  =  2 
irir  =  3 
mj  =  5 


85 

m2  =  5 

ndexing  for  the  stars  on  plate  1  is  f^i^^u 

Star  1  :  the  first  star  on  plate  I,  /iijl  =  lil  =  1.1 

Star  2  :  the  second  star  on  plate  2,  f.ii^l  =  'l^l  ~  2, 1        or      /<r,  ~  - 

Star  3  :  the  third  star  on  plate  1,  fii^l  =  33I  =  3, 1 

Star  4  :  the  forth  star  on  plate  1,  fll^l  =  44I  =  4, 1      or     n^^  =  4 

Star  5  :  the  fifth  star  on  plate  1,  m^l  =  5,5 1  =  5, 1 

ndexing  for  the  stars  on  plate  2  is 

Star  1:  the  first  star  on  plate  2,  y^2i2  =  li2  =  1,2 

Star  3  :  the  second  star  on  the  second  plate,  i.i2j2  =  392  =  3,2 

Star  4  :  the  third  star  on  the  second  plate,  ^232  =  4,32  =  4.2     f.ir2  =  4 

Star  5  :  the  forth  star  on  the  second  plate,  ii2i'2  =  042  =  5,2 

Star  6  :  the  fifth  star  on  the  second  plate^252  =  652  =  6.2     or     /.ij^  =  6 

We  have  four  prime  goals: 

1.  Set  up  the  equations  of  condition. 

2.  Use  these  equations  to  find  J'. 

3.  Use  J',  to  determine  a. 

4.  Use  a,  to  find  /3. 
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Set  up  the  equations  of  condition 

The  equations  F  have  the  form 


F,=  [y]-sJ^"'  )-E,,a,  =  0  (138) 


and  the  equations  G  become 


^       I  ((^cu.    ~  ctu-  )  —cos  6,,^    , 

G=  \c        "  .      X      '  '^      =0,  (139) 


where  k  goes  from  1  to  nir,  the  total  number  of  reference  stars.  The  index  a^^^ 
is  the  /z'*^  number  of  the  k"^  reference  star.  Therefore  in  our  example: 

•  ^n    =  2 

•  ^Ir,     =    4 

It  is  important  to  fix  the  order  of  the  equations.  In  the  last  section  we  stated 
that  the  equations  were  ordered  plate  by  plate  and  within  a  plate,  in  order  of 
increasing  star  number.  However,  from  a  programming  point  of  view  we  found 
it  simpler  to  order  the  equations  star  by  star.  In  other  words  we  group  all  the 
stars  together  that  occur  on  the  separate  plates.  This  is  equivalent  to  ordering  the 
equations  plate  by  plate  (the  row  and  columns  of  the  matrices  are  interchanged). 
Therefore  the  equations  are  ordered  star  by  star  in  order  of  increasing  star  number. 


The  total  vector  H  becomes 
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H- 


-S„a,  \ 


^51 
^51 


Xro 


1^12^2 


iHoi  Si 


'-':?i^i 


>— '32^2 


—  CfiHi 


•—42^ 


Cjiai 


■^5232 


—  £.6232 


y62  -  sr/62 

(0^02  ~  0^2)  ~  COS  ^2 

(^C2    -   ^2) 

(ac4  -  Qi)  -  cos ^4 

(^c6  —  ^e)  —  cos  ^6 
(^c6  -  ^e) 


di2 
621 
CI31 

^32 

d4i 
d42 

d62 
0 
0 

Vo  / 


(140) 


The  (total)  covariance  matrix  is  cr  = 


0 


0 


where  the  covariance  matrix  of 


the  measurements  is 
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I'^ll 

0 

0 

0 

0 

0 

0 

0 

0 

0   \ 

0 

<7l2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

cr^i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(731 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(73, 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

'^41 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

<T^^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

•^51 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

f6i 

0 

V  0 

0 

0 

0 

0 

0 

0 

0 

0 

c^ea/ 

(141) 


Each  sub-matrix  in  cr^  is  a  2  x  2  diagonal  matrix  of  the  form  o-„^  =  (  ^'^  ) . 

n  n 

The  matrix  a^  then  has  dimensions  Y,  -^v  ^  I]  '^^v  The  covariance  matrix 
of  the  reference  stars  has  the  form 


<Ta  = 


(72 

0 

0 

0 

<^4 

0 

0 

0 

c^e 

(142) 


where  the  diagonal  elements  are  2  x  2  block  matrices  of  the  form 


'a/i 


Cifi  COS-  Sfj,        0 


(143) 


so  the  dimension  of  0-3  is  Irrir  x  Inir. 


Bj/  has  dimension  Imu  x  2m,  so  in  this  example  the  total  matrix  B  is  a 
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20x  12  matrix. 


B 


/B„ 

0 

0 

0 

0 

0    \ 

B,. 

0 

0 

0 

0 

0 

0 

B.. 

0 

0 

0 

0 

0 

0 

B31 

0 

0 

0 

0 

0 

B3. 

0 

0 

0 

0 

0 

0 

B„ 

0 

0 

0 

0 

0 

Bj2 

0 

0 

0 

0 

0 

0 

B5> 

0 

0 

0 

0 

0 

B,, 

0 

\  0 

0 

0 

0 

0 

Be  J 

(144) 


The  matrix  K  has  the  dimension  2mrX2m,  in  this  case  it  is  a  6x  12  matrix: 

(145) 


'0     L     0     0     0     0 

K=|0     0     0l2     0     0 

,0     0     0     0     0     L 


Sj.  is  of  dimension  2mi/xl,y  .  Let's  take  a  6-constant  model,  thus  2^  and 
Sj  have  the  dimension  10x6.  Thus  the  total  matrix  2  is  a  20  x  12.  For  each 
star  we  have  a  model  matrix  of  the  form: 


-     ^  f^^.     ^^r     1      0        0      0 
"""       V  0        0       0     ^^,     v^,     1 

and  for  all  the  stars  the  total  model  matrix  is, 


[II 

0    \ 

0 

'—'12 

"21 

0 

■^31 

0 

0 

".32 

241 

0 

0 

—42 

"^51 

0 

0 

■^52 

\  0 

262/ 

(146) 


(147) 


Altogether  the  matrix  A  is 
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A  = 

fsB 

2\ 
0; 

or  expanded 

/sB,, 

0 

0 

0 

0 

0 

"^U 

0    \ 

sBi2 

0 

0 

0 

0 

0 

0 

"12 

0 

sB^i 

0 

0 

0 

0 

•^21 

0 

0 

0 

5B31 

0 

0 

0 

—31 

0 

0 

0 

.SB32 

0 

0 

0 

0 

—32 

0 

0 

0 

.sBji 

0 

0 

"41 

0 

A  = 

0 

0 

0 

,SB,2 

0 

0 

0 

—42 

0 

0 

0 

0 

sBji 

0 

"51 

0 

0 

0 

0 

0 

^652 

0 

0 

"52 

0 

0 

0 

0 

0 

.sB62 

0 

^62 

0 

h 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I. 

0 

0 

0 

0 

\    0 

0 

0 

0 

0 

I2 

0 

0   / 

Deterrr 

line  J' 

(148) 


(149) 


Now  let's  go  back  to  L  =  K^a'^K  +  s~B^a~^B  and  look  at  the  sec- 
ond term:  B^cr~^B.  Since  the  covariance  matrix  is  diagonal,  each  p  di- 
agonal element  is  multiplied  by  the  p"^  column.     The  matrix  dimension  is 

^n2x^0)'^rr'0x '0)^(20x1'^)  ~  ^^^  ^  ^^^"  This  matrix  has  to  compact  blocks  wher- 
ever the  /i'"^  star  was  measured  on  the  v^^  ,  plate. 

/BT<7-ib„+BT<7-ib,2\ 

B2l'721    B21 


B^cr-^B  =  diag 


B 


B3i<73iB3i  +  B'LajB;^^ 
BjiCr^i^Bji  +  B^^cr^j^Bjj 

B62<^6"2   ^62 


42 


V 


(150) 


/ 
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The  matrix  has  dimension  of  K|^j.,^g|(7~g^^  K,(^xi2)  =  (12  x  12) 


K 


T,-l 


K= 


I. 
0 
0 
0 

Vo 

/O       0 

0     a-1 


0 
0 
0 

I. 

0 
0 


0 
0 
0 

vo 


0 
0 
0 
0 


0 
0 
0 
0 

I J 

0 
0 
0 

0     cj 


0 
0 


0 
0 

0 
0 

0 

-1 

! 

0 
0 


a 


0 

,— J 
4 

0 


0 
0 
0 
0 
0 


0    \ 
0 

0 

0 

0 


0    -.-V 


0  I,  0  0  0  0 
0  0  0  I.  0  0 
0     0     0     0     0     I2 


(151) 


Adding  the  two  equations  above  yields  L: 


iT^-li 


|T^-Ii 


L  =  diag 


BuCTjj   Bii   +  Bj2(Tj2  B12  \ 

BT(T-1B,, +a-i 
bI.^.-.^B,,+bT<7-ib,,+^; 


,-1 


\ 


B62'^6~2  B62  +  cr^ 


/ 


(152) 


Using  this  form  of  L  we  can  now  find  J':  J'  =  cr^^  -  s-cr~^BL~^B^cr 
In  our  example,  J'  has  a  dimension  of 


^  2m^  X  ^  2m^  =  20  x  20 


v=\ 


u=\ 


(153) 


From  this  equation  we  see  that  we  first  need  to  find  the  inverse  of  L.  Since  L  is 
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block  diagonal,  its  inverse  is  just  the  inverse  of  the  individual  block  components. 


diag 


/      (BTar/B„+BTa-iB,,)-i 


+  cr: 


v 


So  we  get  BL-^B^ 


B  1  3  L  "'  B  7, 


l].a-/B,. +B>-^B,, 
(BT<T-/B,,+BTa-iB.,,) 


-1 


o  o 


\ 


1\-1 


(154) 


',.      B5.I.; 


'^'■-r^'BL 


BejL^g'Bs',  ' 


and  substituting 

this 

into 

J'  = 

=  '^x 

iBL 

-^B 

w 

yiel 

J'  = 

/Jill 

J  1  13 

0 

0 

0 

0 

0 

0 

0 

0   \ 

J  131 

J  133 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

J  311 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

J  311 

J313 

0 

0 

0 

0 

0 

0 

0 

0 

J331 

J333 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

J  41 1 

J  413 

0 

0 

0 

0 

0 

0 

0 

0 

J  431 

J  433 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Jsii 

J  il3 

0 

0 

0 

0 

0 

0 

0 

0 

J531 

J  533 

0 

I    0 

0 

0 

0 

0 

0 

0 

0 

0 

J633/ 

(156) 


The  notation  above  is  as  follows.  Jaf,c,  a  is  the  fi  number  of  the  star,  b  and  c  are 
plate  numbers.  One  can  think  of  it  as  the  b'''  plate  crossed  with  the  c'''  plate.  For 
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example,  J112  is  the  /(=1  star  on  the  first  plated  crossed  with  the  second  plate. 
As  one  can  see  J212  does  not  exist,  because  the  second  star  is  not  found  on  the 
second  plate.  In  this  example,  this  star  is  a  reference  star.  In  the  equation  above, 
the  main  diagonal  terms 

Kbc  =  ^iL  -  ^'^7ab ^abKa  ^ac(^xac  ( 1 57) 

and  the  off  diagonal  terms  (where  b  does  not  equal  c)  so  ajj^  =  0,  equation 
(157)  becomes 

Jafcc  = -s"'7~a(,Ba6L7a  B„fa-~„^.  (158) 

As  we  stated  earlier,  J'  is  not  block-diagonal.  J'  will  have  a  2  x  2  block  of 
nonzero  elements  only  at  the  intersection  of  those  line  pairs  and  column  pairs 
which  belong  to  stars  that  were  measured  on  the  z/'*'  as  well  as  on  the  /i"^  frame. 

Determine  a 

Now  we  are  ready  to  calculate  the  vector  a.  From  the  overlap  theory  we  have 

2T(j')(d-Sa)  =  0 

sTj'd  =  S^J'Sa  (159) 


S'^J 


■=T    ■[  sT    ■•  -T    T  -jT    T  -T    T  -jT    ^  -T   J  -T    t  -T  j  q 

■S'T    T  ^T    T  n  ■='T    T  sT    1  WT    T  WT   T  =T    t  ^T    t  -jT    t 


(160) 


(162) 
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where 

IT        —■S'Tt         —        4.  ■=''''    J         —  — T   ,         —  _.T    ,         —  —T    ■[         — 

n—  ■S'Tj         •=•        4.'='Tt         —  i^Tj         —  iC-Tt         s 

21   —  ^-i2»^  121^-11   T^  *^3a"'32i»-'3i   T^  •— 43"43i^^4l     I    ^-53«'5ai'^5i 

TT      —  ='''"  T       ■=•      -I-  ~^  f       ~      4.  •=•'■"  1       ■=•      4.  ■=•''"  1       ■=■ 

±ll3    w-iiJ  113^-13    -1-    ^3jJ  313^—33    -r   yw^^*l4i3ta-43    -f-   ^^^^^513^^53 

n       —  "^   1        =■       J.  o"^   T        ■=•       _L  ■=■''''   T        "       -L  '=■'''   T        ??       J.  •=•''"   T        ■=• 

■^■^23    ^-13"133'--13    T^    •—33"  333^-33    "T   '--43J  433^-43    "T   ^-53"533'--53    +   ^-63^633^-63 

S^Jd  is  found  by  multiplying  EJJ  with  the  vector  d. 


=  7,J...d.,=7,    +  Ji.:.<'.=    +=J,J»ii'l3.    +SJ,J3..''3,    +Sj,j3.,d3,   +=J',  Jj,,d,,    +=T   J^.jd^,   +=J,J5,,d5,    +=T   Jj.^d,,  > 

am 

Finally  we  are  able  to  find  the  parameters  a  from 

a=(3^j'Ey^E^j'd.  (164) 


From  our  earlier  discussion  2^J2  does  not  have  the  same  form  as  in  the  single 
plate  case  (i.e.  block  diagonal).  In  this  case  the  matrix  is  banded  bordered  and 
there  are  simplifying  routines  to  invert  this  matrix. 

Determine  3 

Now  let's  find  f3  from 

/?  =  5L-^BT(7-i(d  -  2a)  (165) 
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/■-Il  "ii'ii 


Finally, 


3  = 


m'bL 

" 

-1 

12 

0 

0 

0 

0 

0 

Q 

0 

L 

-1 
22 

''2l''21 

0 

0 

0 

0 

0 

0 

0 

1  -1 uT    .-> 
^33  "31'^3X 

r  -1  t>T    ..-1 
1-33  ''32''32 

0 

0 

0 

0 

0 

0 

0 

Li->7> 

'^41 

•-,-4' 

bT, 

-1 

0 

o 

0 

0 

0 

0 

D 

■-5-5' 

B?,-5-,' 

^55    "52'^52 


d-Sa  = 


di2 

d42 

d^2 


d-Ha  = 


/s„    0  ^ 

0    S12 

S21    0 

=31       0 

0     S32 

/a, 

241    0 

0     S42 
0     2,2 

V    0      =62  y 

Va2 

/dii  -  SiiHi  \ 

Ui2    —   1^12^2 

021  ~  >^2iai 

Uji     —    ligjai 

"32  ~  —32^2 

Q41  ~  "41  ai 

Q42  ~  '^42a2 

"51  ~  '^51^1 

d,2 

—52^2 

\d62  —  262  a2  / 


•^L22B2iCr2i  (d2i  —  £.2iai) 
-'^3^1i^3^^i^3^  -  23181 )  +  ^L3jBj,(7-l(d32  - 

5L-IBT  a-i(d,i  -  24iai)  +  sL-^Bla-.^d,,  - 


-'^55615^51(^51  -  25181)  +  ^Ljs  B^^a^^'id 

^^66  B62<'"6"2  ("^62  -  26282) 


>— 1282 )   \ 

S3282  j 
1—4282  ) 

"5282) 

/ 


(167) 


(168) 


This  concludes  the  example. 


CHAPTER  5 
OBSERVATIONS 


Observations  for  the  Southern  Polar  Zone 

The  Southern  Polar  Zone  plates  were  exposed  between  September  15,  1955 
and  August  7,  1956  at  Sydney,  Australia,  where  the  Yale  southern  astrograph  had 
been  temporarily  mounted.  The  focal  length  of  the  telescope  is  206  cm  and  the 
lens  is  an  8-in  Ross,  stopped  to  4-inches  aperture  to  assure  good  definition  at 
the  plate  edges.  The  measured  plate  area  is  ll°x  11°  on  17-inch  plates  with 
a  plate  scale  of  lOO'Vmm.  An  objective  grating  gave  short  first-order  spectral 
images  three  magnitudes  fainter  than  the  primary,  to  be  measured  for  the  brighter 
stars.  Exposure  times  ranged  from  three  to  ten  minutes.  In  all,  93  plates  were 
exposed  on  the  64  regions,  but  only  the  best  plate  on  each  region  was  selected 
for  measurement.  Dirk  Brouwer  exposed  60  of  the  plates  on  the  4 1  regions.  The 
remainder  of  the  plates  were  exposed  by  Harley  Wood  (2 1  plates)  and  his  staff,  W. 
H.  Robertson  (5  plates)  and  K.P.  Sims  (7  plates).  Brouwer  planned  the  observing 
program  so  as  to  provide  that  in  each  declination  zone  successive  plates  overlap 
50%  in  right  ascension  and  each  declination  zone  also  overlaps  50%  (Hoffleit, 
1971).  This  overlap  pattern  is  illustrated  in  Figure  1.  As  one  can  see  from  the 
figures  this  arrangement  provides  that  a  star  should  occur  on  two  plates  for  the 
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stars  located  at  -70°,  whereas  closer  to  the  pole  a  star  should  be  found  on  a 
minium  of  4  plates.  Some  stars  occur  on  up  to  20  different  plates.  This  strong 
overlap  feature  is  one  the  prime  reason  we  are  reducing  this  catalog. 


Measurement  of  Stars 

The  plates  were  measured  between  1963  and  1965  with  two  different  measur- 
ing machines.  Between  August  1963  and  early  1965,  eighteen  of  the  twenty  plates 
centered  at  -80°  were  measured  with  a  single  screw  Gaertner  engine  (equipped 
with  a  digitizer  and  automatic  IBM  key-punch).  This  machine  was  used  for  most 
of  the  Yale  zone  catalogs  since  1932  and  was  designed  by  Schlesinger  and  de- 
scribed by  him  in  Schlesinger  and  Barney  (1933),  and  the  remaining  46  plates 
were  measured  with  a  2-screw  Mann  machine  installed  in  1965.  This  was  a 
Type  880  Comparator,  built  to  hold  plates  up  to  18  x  18  inches,  equipped  with 
a  DATA  LOGGER  Model  1045  for  automatically  displaying  the  measured  co- 
ordinates to  the  nearest  micron  and  automatically  transferring  the  coordinates  to 
punched  cards  (Hoffleit,  1971). 

Table  6  gives  the  plate  center,  epoch  of  observation,  and  measuring  machine 
used  for  the  64  Yale  plates. 
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Table  4:  The  Yale  Plates 


Plate 

RA 

Dec 

Epoch 
1955.79 

Engine 

Plate 

RA 

Dec 

Epoch 

Engine 

1 

0 

0 

-75 

33 

9 

36 

-80 

1956.05 

G 

2 

1 

0 

-75 

1955.86 

M 

34 

10 

48 

-80 

1956.28 

G 

3 

2 

0 

-75 

1955.93 

M 

35 

12 

0 

-80 

1956.27 

G 

4 

3 

0 

-75 

1955.88 

M 

36 

13 

12 

-80 

1956.28 

G 

5 

4 

0 

-75 

1956.01 

M 

37 

14 

24 

-80 

1956.50 

G 

6 

5 

0 

-75 

1956.03 

M 

38 

15 

36 

-80 

1956.51 

G 

7 

6 

0 

-75 

1956.03 

M 

39 

16 

48 

-80 

1956.30 

G 

8 

7 

0 

-75 

1956.02 

M 

40 

18 

0 

-80 

1956.30 

G 

9 

8 

0 

-75 

1956.05 

M 

41 

19 

12 

-80 

1956.43 

G 

10 

9 

0 

-75 

1956.05 

M 

42 

20 

24 

-80 

1956.78 

G 

11 

10 

0 

-75 

1956.10 

M 

43 

21 

36 

-80 

1956.79 

G 

12 

11 

0 

-75 

1956.28 

M 

44 

22 

48 

-80 

1956.79 

G 

13 

12 

0 

-75 

1956.28 

M 

45 

0 

0 

-85 

1955.79 
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Not  all  star  images  on  the  plates  were  measured.  A  list  of  target  stars  was 
chosen  beforehand.  These  target  stars  to  be  measured  were  chosen  from  the 
Cape  Photographic  Durchmusterung  (CPD,  Gill  and  Kapteyn,  1899)  to  cover 
nine  to  ten  stars  per  square  degree.  Following  a  suggestion  of  Heinrich  Eichhom, 
approximately  3000  mainly  fainter  stars,  chosen  from  the  Melbourne  Astrographic 
Catalog,  were  added  to  the  selection  list  in  order  to  improve  the  astrometric 
accuracy  by  means  of  a  more  through  investigation  of  the  magnitude  effect.  For 
the  stars  north  of  -80°,  including  stars  as  far  north  as  -69°,  all  of  the  stars  in  the 
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Cape  Photographic  Catalogue  for  1950.0,  zones  -64°  to  -80°  (Stoy,  1966)  were 
included.  For  stars  south  of -80°  all  the  stars  in  the  CPD  were  included,  amounting 
to  50%  more  stars  than  are  included  in  the  corresponding  Cape  Photographic 
Catalogue  for  1950.0.  zones  -80°  to  -90°  (Stoy,  1968).  Only  about  10  percent  of 
the  selected  faint  stars  were  bright  enough  to  also  be  in  the  CPD.  According  to 
the  program  originally  envisioned  by  Brouwer,  1 1600  stars  were  to  be  measured. 
The  original  catalog  lists  18702  stars  of  which  2891  are  selected  faint  Melbourne 
stars.  The  limiting  magnitude  of  the  catalog  is  13.6 

Before  the  selected  stars  can  be  measured,  predicted  rectangular  coordinates 
(x,  y)  for  each  star  to  be  measured  were  computed  from  the  right  ascensions 
and  declinations  from  the  source  catalogs  and  the  positions  of  the  plate  centers. 
These  positions  were  precessed  to  the  equinox  1950.0  and  whenever  possible, 
proper  motion  corrections  are  applied  to  bring  the  predicted  positions  into  the  best 
possible  agreement  with  the  observed  coordinates  on  the  plates.  For  half  of  the 
stars  south  of -80°,  no  proper  motions  are  available.  In  case  of  the  stars  selected 
for  the  Melbourne  Astrographic  Catalogues,  no  right  ascension  and  declination 
were  given,  only  the  measured  rectangular  coordinates  on  the  Melbourne  plates, 
and  the  plate  constants.  For  these  stars  the  plate  constants  and  plate  scale  were 
used  to  find  the  rectangular  coordinates. 

On  both  measuring  machines  the  operator  using  the  predicted  coordinates 
searched  for  the  target  stars  and  bisected  the  image,  usually  by  manually  operating 
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the  screws.  Once  a  setting  had  been  achieved,  the  operator  depressed  a  lever 
which  automatically  activated  the  keypunch  to  record  the  screw  reading.  The 
Gaertner  machine  has  a  single  screw  —  so  measures  in  x  and  y  had  to  be  made 
independently,  with  the  plate  rotated  90°  in  the  machine  for  the  y-reading  after 
all  the  x-measurments  were  completed  (Hoffleit,  1967).  The  Mann  machine  is 
a  2  screw  measuring  machine  so  both  coordinates  were  recorded  at  the  same 
time.  Each  plate  was  measured  only  once  in  each  coordinate.  In  accordance 
with  Schlesinger's  precepts  (Schlesinger  and  Barney,  1933)  alternate  regions  (odd 
numbered)  were  measured  "direct",  i.e.,  in  order  of  increasing  right  ascension; 
the  other  (even  numbered)  were  measured  "reversed".  This  was  done  to  reduce 
systematic  errors. 

Many  of  the  selected  faint  stars  were  evidently  too  faint  to  be  measurable 
on  the  Yale  plates.  Occasionally  two  stars  were  found  close  to  the  predicted 
coordinates  for  one  star.  Many  of  these  were  faint  stars.  In  these  cases  both 
positions  were  measured  and  given  the  same  CPD  number  with  an  additional 
designation  of  1  or  2.  Many  program  stars  had  to  be  omitted  because  of 
overlapping  images.  Nearly  100  secondaries  of  double  stars  were  measured  in 
order  to  avert  the  possibility  of  misidentification  of  the  components.  These  stars 
were  also  given  the  same  CPD  number  with  a  designation  of  1  or  2.  Since  the 
plates  were  taken  through  an  objective  grating,  the  bright  stars  in  the  field  all  have 
first-order  spectral  images  that  are  three  magnitudes  fainter  than  the  central  image. 
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In  these  cases  all  3  positions  were  given  the  same  CPD  number  with  a  designation 
of  0  for  the  central  image,  1  for  the  2  primary  images  and  2  for  the  first  order 
images.  Higher-order  spectral  images,  even  when  present,  were  not  measured. 
For  the  fainter  stars  on  our  program,  the  spectral  images  are  not  visible.  The 
measurer  decided  in  every  case  in  the  course  of  measuring  which  spectral  images 
are  suitable  to  be  measured.  From  an  analysis  of  the  percentages  of  stars  of 
successively  fainter  magnitude  for  which  spectral  images  have  been  successfully 
measured  at  the  Yale  engine,  Dorrit  Hoffleit  found  only  a  small  number  of  stars  as 
bright  as  S.O'"  for  which  the  spectra  were  not  measurable,  whereas  very  few  stars 
as  faint  as  9.0"^  have  measurable  spectra.  A  practical  magnitude  limit  seemed  to 
be  8.3*"  (Hoffleit.   1962). 

Magnitude  System 

No  photometry  was  performed  on  the  images,  so  the  magnitudes  of  the  stars 
were  taken  directly  from  the  input  catalogs.  In  fact,  the  magnitudes  as  well  as 
color  indices  and  spectral  types  were  taken  from  the  input  catalog,  the  CPD  and 
the  Cape  Zone  Catalogues.  In  the  case  of  the  faint  stars  where  the  input  catalog 
was  the  Melbourne  catalog  in  which  only  image  diameters  are  given,  the  stars' s 
magnitudes  were  derived  from  the  published  image  diameters  on  the  Melbourne 
plates  and  calibrated  according  to  the  reduction  tables  included  in  the  Melbourne 
Astrographic  Catalogs,  having  a  limiting  photographic  magnitude  of  13.0 


CHAPTER  6 
PLATE  REDUCTIONS 


Reference  Catalog 

As  we  stated  in  Chapter  1,  in  order  to  reduce  photographic  plates  a  set  of 
reference  stars  must  be  used  whose  images  are  on  the  plate.  The  adopted  reference 
catalog  is  the  United  States  Naval  Observatory's  International  Reference  Star 
Catalog  (IRS.  Corbin.  1991).  provided  in  two  versions  FK4/B  1950  and  FK5/J200 
.  The  IRS  is  an  all-sky  catalog  of  positions  and  proper  motions  that  is  based  on  the 
AGK3R  (Corbin  1977,  1978)  in  the  northern  hemisphere  and  the  newly  completed 
SRS  (Southern  Reference  Stars,  Smith  et  al.  1990)  in  the  south.  The  IRS  has 
a  density  of  1  star/deg^  with  most  of  the  stars  primarily  in  the  magnitude  range 
7'".0  to  9'".0.  It  was  compiled  by  matching  122  meridian  circle  catalogs  with 
the  AGK3R  and  the  SRS  to  provide  the  base.  South  of  -60°,  the  IRS,  has  a 
mean  error  in  position  of  0".09  and  0."10  for  right  ascension  and  declination, 
respectively.  The  mean  errors  in  proper  motion  are  ±0.50  arcsec  per  century  and 
±0.5 1  arcsec  per  century  in  right  ascension  and  declination,  respectively.  We  have 
the  fortunate  occurrence  that  the  mean  epoch  of  the  reference  catalog  (  1956.6  for 
right  ascension  and  1954.4  for  declination)  falls  extremely  close  to  the  average 
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epoch  of  the  Yale  observations,  1956.13.  Because  of  this  short  time  interval,  the 
errors  on  the  IRS  proper  motions  can  be  neglected. 

In  order  to  determine  which  images  on  the  plate  are  IRS  stars,  we  applied 
proper  motions  to  the  IRS  stars  to  the  average  epoch  of  Yale  Observations,  deter- 
mined the  standard  coordinates  by  the  gnomonic  projection,  and  then  determined 
preliminary  x  and  y  values  by  using  plate  parameters  determined  from  the  original 
reduction  by  Lu.  However,  while  the  original  solution  used  a  20-constant  model 
(see  equation  175),  we  only  used  the  constants  corresponding  to  the  6-constant 
model.  After  a  preliminary  match  of  the  reference  stars  using  these  approxi- 
mate plate  parameters,  we  performed  a  single  plate  adjustment  determining  more 
accurate  plate  parameters,  and  with  these  improved  plate  parameters  we  then  re- 
matched  the  reference  stars.  A  total  of  1121  IRS  stars  were  found  on  the  plates, 
with  an  average  of  approximately  70  IRS  stars  per  plate.  If  we  were  only  do- 
ing a  single  plate  reduction,  we  would  first  apply  the  proper  motions  to  the  IRS 
stars  to  bring  them  to  the  epoch  of  each  plate,  find  the  matches,  and  perform  the 
adjustment  at  the  epoch  of  each  plate.  However,  this  leaves  the  positions  deter- 
mined from  the  reduction  process  dependent  on  which  plate  the  are  photographed 
on.  A  star  found  on  different  plates  will  have  different  positions  depending  on 
the  epoch  of  the  observation,  again  causing  a  problem  for  the  overlap  solution. 
Since  we  have  a  short  time  interval  between  the  epochs  of  observation  of  Yale 
plates  (maximum  0.83  years)  we  assumed  the  proper  motion  over  this  time  to  be 
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negligible  and  simply  apply  the  proper  motion  of  the  stars  to  bring  them  to  the 
average  epoch  of  the  plates. 

Though  there  are  other  denser  reference  catalogs  — ACRS  and  the  PPM — 
both  of  these  catalogs  are  compilation  catalogs  and  contain  the  original  positions 
from  Yale  catalog,  so  their  use  as  a  reference  catalog  would  have  been  circular 
and  introduced  correlations. 

External  Catalog 

In  order  independently  to  analyze  the  results  we  use  an  external  reference 
catalog,  the  90,000  Star  Supplement  to  the  Positions  and  Proper  Motions  Supple- 
ment Catalog  (we  refer  to  this  catalog  as,  PPM-S,  Roeser,  et  al.,  1993).  These 
stars  are  not  used  in  the  plate  reduction.  We  chose  the  PPM-S  because  it  is 
one  of  the  few  reference  catalogs  in  the  southern  celestial  hemisphere  which  is 
dense  and  accurate  enough  to  give  us  a  meaningful  estimate  of  the  accuracy  of 
our  results  and  which  was  not  complied  using  the  original  Southern  Polar  Zone 
Catalogue.  The  mean  error  in  position  (at  the  average  epoch  of  the  Yale  plates) 
for  those  PPM-S  stars  located  on  plates  is  0."097  in  right  ascension  and  0."094 
in  declination.  The  mean  epoch  for  these  stars  is  1951.64  for  right  ascension  and 
1951.65  for  declination.  Once  again,  we  have  the  favorable  circumstance  that  the 
mean  epoch  of  the  catalog  is  close  to  the  mean  epoch  of  the  Yale  plates,  so  proper 
motion  errors  (0."34  and  0".33,  for  right  ascension  and  declination,  respectively) 
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can  be  considered  negligible.  Using  the  same  procedures  for  the  reference  stars 
we  found  8,768  PPM-S  stars  on  the  Yale  plates,  with  an  average  of  150  PPM-S 
stars  per  plate. 
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Figure  10:  Distribution  of  Magnitudes 


PPM-S  stars  are  used  as  an  external  check  on  the  accuracy  of  our  solution  in 
both  the  single  plate  and  overlap  reductions.  After  a  plate  reduction  we  calculate 
the  right  ascension  and  declination  of  these  stars  using  the  measured  x  and  y 
values,  the  plate  parameters,  and  a  gnomonic  projection.  We  then  compare  these 
calculated  positions  to  the  cataloged  positions,  subtracting  in  quadrature  the  mean 


106 
error  of  the  PPM-S  star  positions.  This  gives  a  independent  check  on  our  solution. 

One  problem  of  using  an  external  catalog  as  a  check  on  the  solution  is  that 
since  the  reference  catalog  is  the  IRS,  the  plate  parameters  establish  a  system  based 
on  the  IRS  reference  system.  There  will  inevitably  be  small  systematic  effects 
between  the  PPM-S  and  the  IRS  reference  systems.  When  using  an  external 
catalog  the  difference  between  the  cataloged  right  ascension  and  declination  and 
the  values  calculated  with  the  plate  constants  will  be  made  up  of  two  parts:  a 
small  part  which  represents  the  difference  between  the  reference  catalog,  and  a 
larger  part  which  determines  how  accurate  our  positions  are.  The  majority  of  the 
stars  in  the  PPM-S  are  fainter  than  those  in  the  IRS.  The  average  magnitude  of  the 
Yale  stars  is  9."'64,  for  the  PPM-S  it  is  Q.'^ST,  and  for  the  IRS  stars  it  is  8.'"3 1 .  In 
fact  the  distribution  of  the  magnitudes  is  very  similar  between  the  PPM-S  and  the 
Yale  stars  (see  Figure  10).  This  is  good  from  the  standpoint  that  most  of  the  Yale 
stars  are  fainter  than  the  IRS  stars  so  we  can  probe  a  magnitude  effect  better  with 
the  PPM-S.  However,  from  the  standpoint  of  comparing  the  results  this  difference 
in  magnitude  may  introduce  a  magnitude  effect  into  our  comparison.  The  PPM-S 
does  not  contain  any  stars  that  are  also  in  the  IRS  catalog,  therefore  it  is  difficult 
to  check  for  systematic  differences  between  the  reference  systems.  Because  the 
PPM-S  is  an  extension  of  the  PPM  system,  we  can  use  the  PPM  (which  contains 
stars  that  also  appear  in  the  IRS  catalog)  as  a  check  of  systematic  differences.  It 
is  probably  safe  to  say  that  the  PPM  and  PPM-S  are  based  on  the  same  system. 
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We  compared  1 804  stars  in  common  between  the  PPM  and  IRS  catalogs  for  the 
-60°  to  -90°  region.  After  proper  moving  both  catalogs  to  the  average  epoch  of 
the  plates  we  found  a  difference  of  -0."0072  and  -0."0049  in  right  ascension  and 
declination,  respectively.  So  we  can  safely  use  the  PPM-S  as  an  external  check 
on  our  solution  since  these  differences  are  at  least  a  factor  of  10  below  the  mean 
error  of  the  PPM-S. 

Determining  the  Plate  Model 

The  choice  of  the  model  is  one  of  the  most  important  steps  in  reducing 
photographic  plates.  With  11°  x  11°  plates  the  relationship  between  the  measured 
and  standard  coordinates  can  be  rather  complex.  The  model  will  invariably  not 
account  for  all  of  the  contributing  effects.  Our  task  is  to  determine  the  model  that 
best  represents  the  physical  situation.  We  start  out  by  making  an  adjustment  using 
a  simple  6-constant  model,  plot  the  residuals  to  check  for  correlations,  and  then 
add  terms  to  the  model  in  an  attempt  to  eliminate  the  correlations.  We  analyzed 
common  optical  effects  such  as  coma  and  optical  distortion,  and  other  possible 
errors  producing  a  non-orthogonal  system,  such  as  plate  tilt  and  rotation.  We 
must  be  careful  when  analyzing  the  terms  to  be  included  in  the  model  and  include 
terms  for  which  we  have  a  "physical"  understanding.  In  Chapter  2  we  pointed 
out  that  main  of  the  effects  we  are  trying  to  model  often  produce  similar  and/or 
correlated  terms  in  the  model  and  therefore  it  is  difficult,  if  not  impossible,  to 
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separate  different  contributions  and  evaluate  the  physical  meaning  of  the  plate 
parameters  in  detail.  In  addition,  we  have  the  problem  of  parameter  variance: 
when  the  number  of  parameters  is  large,  an  unfavorable  error  propagation  of 
the  resulting  star  positions  in  the  plate  field  due  to  the  random  errors  in  the  x,y 
measurements  and  the  reference  star  positions  is  introduced  (Eichhorn  &  Williams 
1963).  Our  approach  for  the  plate  reduction  is  therefore  to  set  up  a  minimal  plate 
model  by  taking  the  physical  meaning  of  the  plate  parameters  into  consideration 
and  applying  known  transformations  and  corrections  to  the  x,y  measurements  a 
priori.  This  minimal  plate  model  then  is  restricted  to  describe  the  properties  of 
the  mapping  and  measuring  process. 

Single  Plate  Reduction 

We  first  determine  the  model  using  single  plate  adjustments.  The  starting 
point  is  the  basic  6-constant  model  (Chapter  2)  and  terms  are  added  to  this  model 
appropriately.  We  perform  a  weighted  least  squares  reduction,  pre-correcting  the 
x,y  positions  of  the  measured  stars  for  refraction.  The  weights  for  the  reference 
stars'  equatorial  coordinates  are  taken  from  the  reference  catalog,  and  for  the 
measurement  errors  we  first  assumed  an  approximation  of  4  microns.  Later  we 
use  a  more  sophisticated  method  for  determining  the  measuring  errors.  Since  the 
single  plate  adjustment  only  uses  the  reference  stars,  reference  stars  with  large 
residuals  (this  could  be  due  to  bad  matches)  will  significantly  effect  the  results. 
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After  the  adjustment  we  determined  the  residuals  using  the  newly  determined  plate 
parameters.  All  reference  stars  with  residuals  greater  than  3  standard  deviations 
from  the  mean  were  eliminated  and  a  new  reduction  was  performed.  In  total 
three  adjustments  were  done.  The  first  two  were  for  the  purpose  of  removing 
high  residual  stars,  and  the  last  reduction  determined  the  final  plate  parameters. 
Using  this  scheme  we  were  able  to  improve  upon  the  precision  of  our  solution 
by  only  removing  a  few  stars.  For  example,  in  the  case  of  the  13-constant  model 
adjustment  we  started  with  4510  measurements,  and  by  the  last  adjustment  we 
had  eliminated  137  measurements,  most  of  these  come  from  a  handful  of  plates, 
and  the  average  number  of  reference  stars  was  68  per  plate.  In  all  iterations,  the 
catalogued  right  ascension  and  declination  values  were  used  in  the  reduction;  we 
did  not  iterate  on  the  solution. 

After  every  reduction  we  used  the  reference  stars  and  the  PPM-S  stars  and 
compared  the  cataloged  right  ascension  and  declination  with  the  computed  values. 
The  average  rms  for  a  single  measurement  for  some  of  the  different  reductions 
we  performed  are  given  in  Table  5.  We  want  to  point  out  that  these  values  are 
actually  the  precision  of  a  single  measurement  and  not  the  averaged  position. 
Since  -  on  the  average-  a  star  is  imaged  on  6.97  different  plates,  we  expect  that 
after  the  average  positions  of  each  star  is  found  these  values  will  decrease  by 
-jr-TT-  So  in  the  case  of  the  13-constant  single  plate  adjustment  the  rms  values 

\/6.9(  •-'1. 

for  average  position  should  roughly  be  (from  the  above  reasoning)  0."22  for  right 
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ascension  and  declination. 

Table  5:  Single  Plate  RMS  Results 
Model  IRS  REF  RMS  PPM-S  EXT 


6-constant 

RA  :  0."  67 

0."  93 

DEC:  0."  66 

l."00 

8-constant 

RA  :  0."  46 

0."  66 

DEC:  0."  49 

0."71 

1 3-constant 

RA  :  0."41 

0."  58 

DEC:  0."  43 

0."  62 

15-constant 

RA  :  0."  40 

0."  55 

DEC:  0."  42 

0."  57 

1 8-constant 

RA  :  0."  39 

0."  59 

DEC:  0."  42 

0."  62 

20-constant 

RA  :  0."  42 

0."  70 

DEC:  0."  44 

0."  72 

Besides  checking  the  cataloged  positions  against  the  computed  positions  we 
determined  the  best  model  by  computing  the  following  residuals  for  the  reference 
stars: 

a;)=(;:)-fe) 

where  x^  and  y^  are  the  measured  x  and  y  values  (pre-corrected  for  refraction) 
and  Xc  and  yc  are  calculated  from  the  cataloged  using  the  plate  parameters.  We  use 
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these  residuals  to  determine  if  we  need  to  add  terms  to  the  model.  The  residuals 
are  plotted  against  possible  terms  which  we  might  want  to  add  to  our  model.  If 
we  find  a  correlation,  then  we  add  the  term  or  terms  to  the  model  and  re-run  the 
plate  reduction,  determining  new  plate  parameters  and  then  replotting  the  graphs. 
The  terms  we  looked  at  are  : 

2        2      /    2    I       2\         /    2    ,       2\         /■    2    ,       2\ 

X,  y,  xy,  X  ,  y  .   (x    +y  ).  x(x    +y  ),  y(x-  +  y-), 

(170) 
m,  mx,  my,  mx",  my",  m",  m"x,  m'y,ci,x'^,y'^,x*y,y"x 

where  m  is  magnitude  and  ci  is  the  color  index.  These  terms  allow  us  to  check 

for  coma,  radial  distortion,  magnitude  effects,  and  other  second  and  third  order 

effects. 

Figures  14-18  are  the  residuals  for  the  6-constant  model  for  all  the  reference 
stars  on  all  the  plates  (after  a  3  sigma  elimination  test).  A  third  degree  polynomial 
was  fit  to  the  points,  for  the  purpose  of  analyzing  the  trends  in  the  data.  The  left- 
hand  side  of  the  plots  are  for  dx  and  right-hand  side  for  dy.  From  these  figures 
we  can  make  the  following  conclusions:  there  is  a  strong  correlation  of  dy  with 
y,  there  is  a  strong  indicator  for  the  inclusion  of  a  plate  tilt  term  because  dx  is 
correlated  with  both  xy  and  x^  and  dy  is  correlated  with  xy  and  y^,  and  there 
also  is  a  good  indication  that  a  distortion  term  should  be  included  because  of 
the  correlation  of  both  dx  and  dy  with  x(x^  -f-  y^),y(x^  +  y~)-  Because  of  the 
complicated  way  the  residuals  interact  we  will  say  nothing  else  about  the  other 
plots  until  we  have  included  the  above  terms. 
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In  Figures  19-23  we  included  the  terms  only  for  plate  tilt.  This  model  is  an 
8-constant  model  has  the  form 

'<f     '?     1     0     0     0     (~     ^T] 


yO    0    0    ^    rj    I    e^/    rj'  ^^^^^ 

This  dramatically  improved  the  results  and  demonstrates  the  cross  interaction  these 
terms  have  on  one  another.  However,  there  is  still  a  radial  distortion  term  that 
should  be  included.  Using  a  13-constant  model  (equation  172)  we  included  the 
radial  distortion  term  and  split  the  second  order  effects  for  both  dx  and  dy,  see 
figures  24-28.  This  model  has  the  form: 
2^/^^  r;  1  0  0  0  ^2  ^^  rj^-  0  0  0  ^{e  +  n')\  ^,^^^ 
VO  0  0  ^  ^  1  0  0  0  e  ^V  r  rjie  +  rj-'))-  ^^^^^ 
Though  we  added  a  term  for  radial  distortion  the  residual  plots  for  this  term  still 

show  a  systematic  trend.  However,  this  is  probably  due  to  the  interaction  of  the 
x^y  and  y^x  terms.  After  adding  these  terms  to  the  13-constant  model  we  have 
a  15-constant  model  of  the  form 

^^[^    V    ^    0   0   0   e    ^v   r    0     0     0    ^ie  +  r)    a^v)      o    \ 

VO     0     0     ^     r?     1      0       0       0      e     ^V     r     rj{e+r)         0         riiir]))' 

(173) 

From  the  residual  plots  (figures  29-33)  we  have  graphed,  we  see  that  all  the 

correlations  are  essentially  eliminated: 

We  tested  many  different  models,  including  an  1 8-constant  model  which 
added  first-order  magnitude  terms  to  the  15-constant  model: 

^£     r,     1     0     0     0     {2     {,,     ,,2      0       0       0      {({2  +  7,2)     {({,,)         0         (m-mo)  0  {(m-mo)> 


0    0    0    i    n    1     0     0      0     e    in    n''    nie  +  n")       o       ,j({,j)  o  (m-mo)    n(m-mo)/ 
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and  the  original  20-constant  model  adopted  by  Lu: 


/«      1     '     0     0     0     £-      Si     I-*       0        0        0       (m  -   mo)  0  (m-mo){  0  (c-C(,){  0  (m-mol.| 

\0     0     0     £      1     I       0        0        0       £^     Si,      ,^  0  (m  -  mo)  0  (m-mo)i)  0  (c-co)"!  0 


(m  -  mo)« 


However,  we  have  only  shown  the  residual  plots  for  20-constant  model  Figures 
34-38).  We  determined  that  the  "best"  model  for  the  single  plate  solutions  was  the 
15-constant  model.  From  Table  5  we  can  see  that  this  model  gives  the  smallest 
rms  and  basically  removes  all  the  systematic  trends  in  the  residual  plots.  However, 
the  difference  between  this  model  and  the  18-constant  model  is  quite  small.  The 
main  difference  between  these  two  models  is  with  the  magnitude  terms.  The 
magnitude  range  of  the  IRS  catalog  is  roughly  2  magnitudes  (T'^-Q'"),  [see  Figure 
10],  and  with  this  small  range  it  is  difficult  to  test  for  a  magnitude  term.  The 
situation  will  improve  for  the  overlap  solution,  where  all  the  stars  are  used  in  the 
reduction,  so  our  magnitude  range  increases  to  5  magnitudes  (7'"-12'")  and  the 
usage  of  the  grating  images  allows  us  to  explore  magnitude  effect  in  greater  detail. 
For  this  reason  and  the  desire  to  have  as  few  plate  parameters  as  needed  we  chose 
the  15-constant  model  to  determine  the  positions  for  the  single  plate  solution. 
We  found  that  our  model  fit  the  data  better  than  the  20-constant  model  adopted 
by  Lii.  We  did  not  find  any  justification  for  the  inclusion  of  the  magnitude  nor 
color  terms. 
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Final  Positions  from  the  Single  Plate  Solution 

As  we  stated  before,  the  major  problem  with  the  single  plate  approach  is  that 
a  star  located  on  n  different  plates  has  n  different  position  estimates.  In  Figure 
1 1  we  demonstrated  this  by  comparing  the  positions  from  plates  1  and  plate  24 
(top)  and  plates  1  and  25  (bottom)  after  a  single  plate  reduction  using  the  plate 
constants  from  a  15-constant  model.  In  Figure  12  we  show  three  of  the  polar 
plates,  plate  61  overlapped  with  plate  64  and  plate  62  overlapped  with  63.  In  the 
bottom  left  hand  comer  a  scale  of  2  arcseconds  is  given.  The  values  given  on 
the  upper  right  are  the  average  standard  deviation  in  the  positions  from  the  two 
different  plates.  RMS  ALL  is  the  average  difference  in  the  position.  In  order 
words,  it  is  the  average  length  of  the  lines.  The  plates  overlapped  with  plate  1 
have  on  the  average  about  0."8  difference  in  position,  while  for  the  polar  plates  we 
have  closer  to  1".  Though  these  differences  appear  large,  these  are  the  individual 
differences  in  position  from  one  plate  to  the  other  and  not  the  difference  between 
the  individual  value  and  the  average  value. 
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Figure  11:  The  different  in  positions  found  from  the  single 
plate  adjustment  for  plates  1  and  24  and  plates  1  and  25. 
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Figure  12:  The  different  in  positions  found  from  the  single 
plate  adjustment  for  plates  61  and  64  and  plates  62  and  63. 
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After  determining  all  the  stars'  equatorial  coordinates  (using  the  plate  pa- 
rameters from  the  single  plate  adjustment)  we  then  matched  the  images  of  the 
same  star  on  different  plates  and  averaged  these  positions  together.  This  step  is 
probably  the  most  frustrating  and  time  consuming.  The  data  for  the  Yale  plates 
were  measured  on  semi-automatic  measuring  machines  and  much  of  the  informa- 
tion concerning  the  observations  and  measuring  has  been  lost.   A  great  deal  of 
time  (and  patience)  was  required  in  cleaning  up  this  data  and  matching  the  star 
images.  When  the  plates  were  measured  there  were  hundreds  of  "possible  target 
stars"  which  were  close  to  the  program  stars  and  therefore  were  measured  in  ad- 
dition to  the  probable  target  star  to  avert  any  confusion.  Part  of  the  "cleaning" 
process  involved  determining  which  of  these  measurements  were  truly  images  of 
the  intended  program  stars.  We  also  had  to  eliminate  other  images  either  because 
they  were  imperfections  on  plate,  bad  measurements  or  repeated  measurements. 
In  addition  there  were  numerous  cases  where  only  one  measurement  was  given 
or  where  the  magnitude  of  a  star  was  definitely  incorrect.  We  deleted  the  images 
with  only  one  measurement  and  corrected  the  magnitudes  by  cross  referencing 
them  with  the  CPD  catalog. 

The  choice  for  the  matching  tolerance,  of  course,  determines  which  images 
we  call  the  same  star.  It  is  extremely  important  to  correctly  match  the  images  of 
the  same  stars  from  different  plates.  If  there  is  too  low  a  matching  tolerance,  then 
not  all  the  images  of  the  same  star  will  be  matched  together  and  this  will  result 
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in  artificially  increasing  the  total  number  of  stars  we  determine.  If  we  have  too 
high  a  tolerance  then  "bad"  data  and  images  which  really  belong  to  other  stars 
will  be  combined  together.  This  will  artificially  decrease  the  total  number  of  stars 
we  determine  and  weaken  the  solution. 

In  total  there  are  124.229  star  images  of  17506  stars.  We  found  3831  stars 
were  averaged  together  because  they  were  measured  twice  by  accident  and  1600 
possible  double  stars  which  were  removed.  We  do  not  want  to  include  double  stars 
in  the  overlap  solution  so  we  removed  them  and  treat  these  stars  separately.  There 
were  a  number  of  images  that  were  on  only  one  plate  (2,752)  either  because  they 
occurred  at  the  borders  of  the  Southern  Polar  Zone  (so  they  are  only  imaged  on 
one  plate)  or  they  were  simply  bad  data  points.  These  measurements  are  not  used 
in  the  overlap  solution.  Their  positions  are  determined  from  the  plate  parameters 
after  the  overlap  solution. 

After  cleaning  the  data  and  matching  and  averaging  the  images  together  we 
compared  the  average  position  with  the  individual  values  from  the  plates  and 
found  the  rms  in  right  ascension  to  be  0."54  and  in  declination  0."55.  These 
values  are  the  precision  of  the  individual  measurements  and  are  consistent  with 
the  values  from  Table  5. 

We  also  compared  our  averaged  single  plate  position  with  those  found  in 
the  PPM-S  catalog.    For  1874  stars  in  common  we  obtained  an  rms  in  right 
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ascension  of  0".28  and  in  declination  0".32.   After  quadratically  subtracting  the 
error  due  to  the  PPM-S  we  get  0."26  and  0."31,  for  right  ascension  and  declination, 
respectively.  These  values  are  roughly  the  mean  error  in  our  positions  after  the 
single  plate  solution. 

Overlapping  Plate  Reduction 

As  we  stated  earlier  an  overlap  adjustment  uses  all  the  stars  that  either  occur 
on  more  than  one  plate  or  are  reference  stars.  This  dramatically  increases  the 
number  of  observations  used  in  the  adjustment.  In  the  case  of  the  single  plate 
method  we  had  roughly  70  reference  stars  per  plate;  however,  now  we  have 
124,229  star  images  of  17,506  stars  greatly  increasing  the  computational  effort. 

Besides  using  the  PPM-S  stars  as  a  means  of  evaluating  the  accuracy  of 
our  solution  we  also  randomly  pulled  approximately  10  IRS  stars  on  each  plate 
which  were  flagged  not  to  be  used  as  a  reference  stars.  So  as  not  to  confuse 
these  stars  with  IRS  reference  stars  we  call  these  flagged  stars  "external"  IRS 
stars.  After  every  reduction  we  use  the  reference  stars,  external  IRS  stars  and  the 
PPM-S  stars  and  compare  the  cataloged  right  ascension  and  declination  with  the 
computed  ones.  The  average  for  all  the  plates  of  the  rms  values  for  some  of  the 
different  reductions  we  performed  is  given  in  Table  6 

As  in  the  single  plate  analysis  we  tested  various  different  models,  starting  with 
the  6-constant  model  and  adding  terms  appropriately.  In  this  case  we  iterated  on 
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the  solution  three  times;  after  every  iteration  we  added  the  correction  to  the 
positions  to  the  field  stars,  but  not  the  reference  stars.    In  addition,  instead  of 
assuming  a  measuring  error  of  4  microns  we  used  a  measuring  error  dependent 
on  magnitude.  The  measuring  error  associated  with  a  measuring  machine  depends 
on  how  precisely  the  measuring  machine  can  determine  the  center  of  the  image. 
If  the  image  is  large  (i.e.,  a  bright  star)  or  if  the  image  is  very  small  (i.e.,  a  faint 
star)  the  center  is  difficult  to  determine.  Therefore,  bright  and  faint  stars  generally 
have  a  larger  measuring  error  associated  with  them.   Every  measuring  machine 
has  its  own  characteristic  range  of  image  diameters  where  its  measuring  error 
is  at  a  minimum.  In  order  to  find  a  relationship  to  characterize  the  magnitude- 
dependent  measuring  error  we  used  the  stars  in  the  PPM-S  and  the  PPM.  Here  we 
have  included  the  stars  in  the  PPM  because  they  are  brighter  than  the  PPM-S  and 
we  want  to  adequately  cover  the  entire  magnitude  range  found  on  the  plates.  We 
calculated  the  predicted  x  and  y  values  for  the  stars  in  these  catalogs  and  compared 
them  to  the  measured  x  and  y  values  and  subtracted  the  error  in  the  positions  from 
the  reference  catalog.  We  have  assumed  that  most  of  the  remaining  error  is  due  to 
measuring  error  and  not  systematic  errors  remaining  from  incorrect  modeling.  In 
this  case,  we  used  the  plate  parameters  that  included  terms  for  linear  magnitude 
effects  so  that  our  results  would  not  be  corrupted  by  a  magnitude  term.   Since 
the  plates  were  measured  on  two  different  measuring  machines  we  performed  the 
analysis  separately  for  each.  We  fit  a  third-degree  polynomial  through  the  points 
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and  used  this  as  our  function  to  estimate  the  measuring  error  for  the  overlap 
reduction.  The  graph  of  these  results  is  shown  in  Figure  13  (after  the  points  were 
binned).  The  binned  points  are  shown  as  circles  —  with  large  circles  indicating 
few  points  are  contained  in  the  final  binned  point  and  small  circles  indicating 
a  large  number  of  points  make  up  the  binned  point.  The  function  representing 
the  magnitude-dependent  measuring  error  for  the  two  machines  is  roughly  the 
same.  However,  for  determining  the  measuring  error  for  the  overlap  we  used  the 
function  corresponding  to  the  appropriate  measuring  machine. 

One  other  change  we  made  for  the  overlap  reduction  was  the  inclusion  of 
the  grating  images  for  the  purpose  of  eliminating  magnitude  effects.  There  are 
6901  stars  with  grating  images.  With  ideal  magnitude-independent  imaging,  the 
primary  stellar  image  will  be  at  the  average  of  the  first  order  images  positions. 
Coma  and  magnitude  effects  will  cause  deviations  from  an  exact  agreement  of  the 
average  of  the  grating  images  and  the  primary  image.  By  examining  this  deviation 
we  can  model  the  magnitude  dependent  terms.  The  stars  with  grating  images  are 
brighter  than  the  average  field  stars  (m<9.0)  and  the  grating  images  themselves 
are  three  magnitudes  fainter  then  the  primary.  Therefore  we  can  use  the  grating 
image  information  to  ensure  that  any  magnitude  effects  are  equally  modeled  for 
the  whole  stellar  magnitude  range  (i.e.  9"^  to  12"^  ).  We  do  this  by  averaging  the 
X  and  y  positions  of  the  grating  images  and  assign  a  fictitious  star  to  the  average 
which  is  also  three  magnitudes  fainter  than  the  primary  and  has  the  same  right 
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ascension  and  declination  as  the  primary.  We  treat  data  for  these  fictitious  stars 
as  additional  information  in  the  overlap,  but  in  the  theoretical  framework  of  the 
overlap  solution  they  are  treated  as  the  same  star  as  the  primary. 
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Figure  13:  Magnitude  Dependent  Measuiing  Error, 
the  units  of  the  variances  are  arcseconds  squared. 


124 
Table  6  lists  the  results  of  some  of  the  models.  The  columns  2— ^  are  the 
rms  values  of  comparing  the  position  determined  from  the  plate  parameters  to  the 
position  from  a  reference  catalog.  The  reference  catalogs  are  as  follows:  IRS 
REF  the  IRS  stars  used  in  the  reduction,  IRS  EXT,  the  IRS  stars  not  used  in 
the  reduction.  PPM-S  IND  the  PPM-S  compared  to  the  individual  value,  and 
PPM-S  AVE  the  PPM-S  compared  to  the  average  position.  The  last  columns  is 
the  comparison  of  a  star  to  its  individual  plate  values. 
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Table  6:  Overlap  Plate  RMS  Results 

Model  IRS  REF  IRS  EXT       PPM-S  IND     PPM-S  AVE    IND  to  AVE 

6-constant       RA:  0."  84  0."  72  0.87"  0.51"  0."72 

DEC:  1.32  1  "  13  1".48  1"22  0."  76 


8-constant       RA:  0."  50  0."  41  0."  61  0."  33  0."  52 

DEC:  0."  54  0."  33  0."  66  0."  37  0."  56 


13-constant      RA:  0."  45  0."  38  0 "  53  0."  28  0."  46 

DEC:  0."48  0."  34  0."  56  0."  32  0."  49 


15-constant      RA:  0."  42  0."  35  0."  50  0."  24  0."  43 

DEC:  0."43  0."  31  0."  53  0."  29  0."  46 


18-constant      RA:  0  "  41  0  "  35  0."  49  0  "  24  0  "  43 

DEC:  0."43  0"  31  0."  51  0 "  27  0."  46 
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Since  in  the  single  plate  adjustment  each  solution  was  dependent  on  only 
60-70  observations,  the  solutions  could  be  somewhat  loose  and  therefore  some 
of  the  systematic  effects  could  be  hidden  in  the  noise.  In  the  case  of  an  overlap 
solution  the  restriction  that  a  star  have  the  same  position  from  plate  to  plate 
makes  the  solution  much  stronger  and  the  systemic  effects  more  noticeable.  The 
residual  plot  for  the  overlap  solution  only  contains  the  weighted  least  squares 
curve  and  not  the  actual  data  points  (with  124,229  images  the  graphs  would  have 
been  black  bands).  We  did  not  include  the  residual  plots  for  the  6  or  8-constant 
model.  The  systematic  trends  for  these  graphs  are  similar  to  the  single  plate 
results:  after  a  reduction  with  a  6-constant  model  the  residual  plots  showed  a 
strong  indication  for  the  inclusion  of  full  second  order  terms  and  radial  distortion. 
The  8-constant  model  improved  the  situation  but  a  distortion  term  still  needed 
to  be  included.  Following  the  procedure  from  the  single  plate  we  tested  the 
main  models  including  our  13-constant  (Figures  39^2)  15  (Figures  43-46)  and 
1 8-constant  models  (Figure  47-50). 

Conclusions 

From  Table  6  we  can  see  that  the  15-constant  and  1 8-constant  model  give 
roughly  the  same  values.  However,  we  feel  that  the  coma  term  and  linear 
magnitude  terms  should  be  included.  First  of  all,  it  is  very  likely  that  with 
such  large  plates  (11°  x   11°)  a  coma  term  ,  though  small,  exists.   In  addition, 
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technically  since  we  have  included  the  grating  images  as  additional  information, 
not  including  a  magnitude  term  like  coma  will  probably  add  to  the  noise  of  the 
other  parameters.  A  linear  magnitude  term  corrects  for  guiding  error  which  can 
vary  greatly  from  plate  to  plate.  Though  this  coefficient  might  be  near  zero  on 
most  of  the  plates,  it  is  wise  to  include  it  in  case  it  is  large  on  a  few  plates.  Though 
we  are  not  as  limited  in  the  single  plate  case  to  keep  the  number  of  parameters 
to  the  absolute  minimum,  including  too  many  terms  will  cause  problems  in  our 
inversion  process,  leading  to  either  an  non-invertible  matrix  or  at  best  an  poorly 
determined  matrix.  From  Table  6  and  the  residual  plots  we  decided  the  "best" 
model  was  the  18-constant  model.  Using  this  model  we  determined  the  positions 
for  all  the  stars  and  in  Tables  7  and  8  we  compared  our  results  for  the  single 
plate  and  overlap  solutions. 

Table  7:  RMS  Single  Plate  Comparisons:    15-constant  model 

RA  (")  DEC(") 

IRS  REF  TO  IND  0.40  0.42 

PPM-S  TO  IND  0.55  0.57 

PPM-S  TO  AVE  0.27  0.29 

IND  TO  AVE  0.54  0.55 
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Table  8: 

RMS 

Overlap 

Plate  Comparisons: 
RA  (") 

18-constant  model 
DEC  (") 

IRS  REF  TO  AVE 

0.41 

0.43 

PPM-S  TO  IND 

0.49 

0.51 

PPM-S  TO  AVE 

0.24 

0.27 

IND   TO  AVE 

0.43 

0.46 

The  18%  improvement  in  the  precision  of  the  average  position  over  that 
the  individual  measurement  between  the  overlap  and  single  solutions  shows  an 
improvement  in  the  precision  of  positions  due  to  a  stronger  adjustment  producing 
better  determined  plate  parameters.  In  the  original  reduction  quoted  this  value  at 
0."55.  The  original  reduction  used  the  Second  Cape  Catalogue  as  the  reference 
catalog.  Using  this  catalog  there  was  approximately  100-200  stars  per  plate.  We 
have  on  the  average  70  stars  per  plate.  Thus,  though  we  have  a  better-determined 
reduction  model  the  results  are  roughly  the  same  between  the  two  reductions 
because  of  the  dense  reference  catalog  used  in  the  original  reduction.  There  is 
an  improvement  of  1 1  %  from  the  single  plate  reduction  and  the  overlap  in  the 
comparison  of  the  PPM-S  to  the  average  vaule.  This  improvement  would  be  larger 
if  the  effect  of  parameter  variance  were  removed  In  addition,  the  measuring  error 
of  4  microns  is  now  dominate  and  it  is  difficult  to  obtain  much  more  improvement. 
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Future  Work 

The  positions  estimates  derived  in  this  investigation  would  be  improved  if 
the  plates  were  remeasured.  We  have  thoroughly  tested  the  improvement  in  the 
precision  attainable  by  the  overlapping  plate  reduction  method.  An  interesting 
addition  to  the  reduction  method  would  be  to  restrict  the  plate  parameters  which 
do  not  vary  much  from  plate  to  plate  to  a  constant  value.  Alternately,  these 
plate  parameters  could  be  constrained  to  an  average  value  which  is  only  allowed 
to  varying  according  known  distribution,  (i.e.  stochastically  constrained).  In 
theory,  by  scholastically  constraining  the  parameters  we  can  find  estimates  for 
parameters  that  would  be  more  accurate  than  those  obtained  by  ignoring  the  fact 
that  they  are  samples  of  a  population  with  known  characteristics.  For  example, 
if  the  plates  were  all  taken  with  the  same  telescope  then  we  would  expect  the 
plate  parameters  corresponding  to  a  scale,  coma,  and  distortion  could  possibly  be 
restricted  to  a  constant  value  or  constrained  to  an  average  value  which  is  allow 
to  vary  according  to  a  known  distribution. 

Though  these  are  interesting  additions  to  the  adjustment  method  and  we  did 
test  constraining  some  of  the  parameters  to  constant  values,  this  data  set  is  simply 
not  accurate  or  precise  enough  to  warrant  such  additions.  However,  if  the  plates 
were  remeasured  or  if  we  used  data  from  one  of  the  current  measuring  machines 
with  a  measuring  error  of  1  micron  then  the  incorporation  of  these  constraints 
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might  lead  to  a  substantial  improvement.  We  considered  remeasuring  the  plates. 
However,  because  of  the  large  size  of  the  plates  and  because  of  the  large  number 
of  them  this  has  proved  beyond  the  scope  of  this  dissertation.  There  are  very  few 
measuring  faculties  that  could  be  used  for  this  project.  Hopefully,  in  the  future 
these  plates  can  be  remeasured. 


The  obvious  next  step  is  to  determine  the  proper  motions  for  the  stars  located 
on  the  plates.  The  best  usage  of  this  newly  compiled  catalog  is  as  a  mid-epoch  or  a 
first  epoch  source.  The  Astrographic  Catalogue  has  been  transferred  to  machine- 
readable  form  and  would  serve  as  an  excellent  first  epoch  with  either  the  CPC- 
2  (the  Second  Cape  Photographic  Catalog,  Zacharias,  1992)  or  the  FOCAT-S 
(Bystrov  et  al.,  1989)  catalog  as  the  third  epoch. 
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Figure  14:  6-Constant  Single  Plate  Residuals 
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Figure  15:  6-Constant  Single  Plate  Residuals 
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Figure  16:  6-Constant  Single  Plate  Residuals 
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Figure  17:   6-Constant  Single  Plate  Residuals 
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Figure  19:   8-Constant  Single  Plate  Residuals 
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Figure  20:  8-Constant  Single  Plate  Residuals 
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Figure  21:   8-Constant  Single  Plate  Residuals 
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Figure  22:  8-Constant  Single  Plate  Residuals 
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Figure  23:  8-Constant  Single  Plate  Residuals 
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Figure  24:   13-Constant  Single  Plate  Residuals 
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Figure  25:   13-Constant  Single  Plate  Residuals 
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Figure  26:   13-Constant  Single  Plate  Residuals 
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Figure  27:   13-Constant  Single  Plate  Residuals 


145 


S 

nqle  Plat 

5  Resid 

UOIS 

13 

Constont  Mode 

20 
10 
0 
10 

1 

1 

i 

■ 

t 

■ 

1   • 

20 

'. 

-0.167  0.277  0.720  1.163 


-20 
-0.167 


0.277  0.720  1.163  1.606 


Figure  28:   13-Constant  Single  Plate  Residuals 
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Figure  29:   15-Constant  Single  Plate  Residuals 
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Figure  30:   15-Constant  Single  Plate  Residuals 


148 


Single  Plote  Residuols  15  Constant  Model 


Wmm^WW^^^^^'': 


-20  [ . 

-2000 


0 
my 


2000  -2000 


0 
my 


1.13«I0  2.25x10''  3.38x10  4.50x10 


■        :M  j    I  i  '  i  ^  .: 


%■   ■•  .y-    }  f.    ■<.-. 


0 
m'x 


20 
10 
0 
10 
20 

■    -^vi^^g^i^^ 

Ba'^ 

:    ~^^^^^^^^ 

!■.'■'-■'  '  •*.  ■  •  ■'           1 

0 
m'x 


Figure  31:   15-Constant  Single  Plate  Residuals 
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Figure  32:   15-Constant  Single  Plate  Residuals 
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Figure  33:   15-Constant  Single  Plate  Residuals 
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Figure  34:  20-Constant  Single  Plate  Residuals 
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Figure  35:  20-Constant  Single  Plate  Residuals 
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Figure  36:  20-Constant  Single  Plate  Residuals 
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Figure  37:  20-Constant  Single  Plate  Residuals 
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Figure  38:  20-Constant  Single  Plate  Residuals 
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Figure  39:   13-Constant  Overlap  Residuals 
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Figure  40:   13-Constant  Overlap  Residuals 
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Figure  41:   13-Constant  Overlap  Residuals 
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Figure  42:   13-Constant  Overlap  Residuals 
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Figure  43:   15-Constant  Overlap  Residuals 
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Figure  44:   15-Constant  Overlap  Residuals 
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Figure  45:   15-Constant  Overlap  Residuals 
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Figure  46:   15-Constant  Overlap  Residuals 
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Figure  47:   18-Constant  Overlap  Residuals 
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Figure  48:    18-Constant  Overlap  Residuals 
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Figure  49:   18-Constant  Overlap  Plate  Residuals 
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has  culminated  with  this  dissertation.  Jane  took  an  extended  hiatus  from  Florida, 
traveling  with  her  family  to  the  Keystone  State.  During  these  formative  years  in 
Pennsylvania  the  Morrison  clan  represented  the  United  States  as  spectators  in  the 
1972  Montreal  Olympic  Games.  After  failing  to  medal  the  dejected  Morrisons 
soon  picked  up  stakes  and  moved  to  Music  City  USA  to  pursue  other  interests.  On 
a  family  excursion  to  the  Great  Smoky  Mountains,  downwind  from  the  Oak  Ridge 
National  Lab,  Jane  was  bitten  by  a  radioactive  spider  which  gave  her  temporary 
prescient  abilities,  allowing  her  to  foresee  her  life  long  ambition  in  mathematics 
and  the  sciences.  The  Tennessee  educational  system  offered  itself  as  a  springboard 
to  accomplish  these  goals.  Disguised  as  a  mild-mannered  geeky  sUident,  voted 
"Most  Studious"  in  Brentwood  High  School,  Jane  graduated  and  went  on  to 
Middle  Tennessee  State  University  where  she  studied  under  the  eminent  Professor 
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Robert  Carlton.  On  an  unusually  blustery  spring  day  in  1987  she  received  her 
Bachelors  of  Science  degree.  Following  the  Sun  south  in  the  Indian  summer  of 
the  fall  of  '88,  Jane  returned  to  the  Sunshine  State  to  fulfill  her  destiny.  After 
thoroughly  exhausting  the  possibilities  in  the  New  World,  Jane  now  sets  her  sails 
for  the  Old  World  to  seek  out  new  adventures. 
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